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Abstract 

IceTop, the surface component of the IceCube Neutrino Observatory at the South Pole, is an air shower array 
with an area of 1 km 2 . The detector allows a detailed exploration of the mass composition of primary cosmic 
rays in the energy range from about 100 TeV to 1 EeV by exploiting the correlation between the shower energy 
measured in IceTop and the energy deposited by muons in the deep ice. In this paper we report on the technical 
design, construction and installation, the trigger and data acquisition systems as well as the software framework 
for calibration, reconstruction and simulation. Finally the first experience from commissioning and operating the 
detector and the performance as an air shower detector will be discussed. 

Keywords: IceCube, IceTop, cosmic rays, air shower, detector 
PACS: 95.55.Vj 
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IceCube Lab 




Figure 1: The IceCube detector with its components DeepCore and IceTop in the final configuration (December 2010). 

1. Introduction 

The IceCube Neutrino Observatory at the geographic South Pole consists of a cubic-kilometer detector situated 
in the ice at a depth between 1450 m and 2450 m and a square-kilometer detector array at the surface. IceCube 
is primarily designed to measure neutrinos from below, thus using the Earth as a filter to discriminate against 
background induced by cosmic rays H] 12 O. The detector employs optical sensors to detect light of charged 
particles generated by neutrinos in the ice or the Earth's crust. In addition, IceCube also includes a more densely 
instrumented part called DeepCore and an extensive air shower array on the surface called IceTop, both fully 
integrated into the IceCube data acquisition system. The IceTop array extends IceCube's capabilities for cosmic 
ray physics allowing the use of the full IceCube Observatory (Fig. [TJ as a 3-dimensional array for the study of 
high-energy cosmic rays. Construction of IceCube, including the IceTop component, was completed in December 
2010. 

IceTop consists of Cherenkov tanks filled with clear ice that operate on the same principle as the water tanks 
of the Haverah Park experiment [4 | and the Pierre Auger Observatory @. The tanks are arranged in pairs on the 
same, approximately 125 m, triangular grid as the vertical cables that carry the deep sensors of IceCube. The two 
tanks at each surface station are separated from each other by 10 m. Each tank contains two standard IceCube 
digital optical modules (DOMs, see Section [3TT] >. Air showers initiated in the atmosphere by cosmic rays are 
typically spread over a number of stations. The light generated in the tanks by the shower particles (electrons, 
photons, muons and hadrons) is a measure of the energy deposit of these particles in the tanks. The information 
from multiple stations is fitted to a model of the overall shower shape and intensity, called 'shower size' S , and the 
direction described by the zenith and azimuth angles 9, <p. 

With the IceTop detector configuration at the 2835 m altitude of the South Pole surface, the threshold for 
efficient reconstruction of air showers that trigger 3 or more stations is approximately 300 TeV in most of the 
detector and approximately 100 TeV in a denser in-fill region. The geometrical acceptance of the combined surface 
and in-ice detectors is approximately 0.3 km 2 sr. This aperture will provide a useful rate of showers with energies 
up to few EeV. Thus, the energy range of IceCube as a cosmic -ray detector fully covers the knee region of the 
spectrum and extends to the energies where a transition from galactic cosmic rays to a population of extra-galactic 
particles may occur. The key to identifying a transition from one population to another is to know the nuclear 
composition of the primary cosmic radiation. General reviews of the status of cosmic ray physics in the region 
covered by IceTop can be found, for example, in (6) or in the cosmic ray review in Q. A recent summary of 
composition measurements is given in |S). 
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Table 1: Acronyms used in this article. 



ADC Analog-to-Digital Converter 

ATWD Analog Transient Waveform Digitizer (front-end chip for signal recording) 

DOM Digital Optical Module (light sensor and readout electronics ) 

DOMHub Computer in the IceCube Lab which receives the signals from the DOMs 

fADC fast Analog-to-Digital Converter (used in the DOM to record waveforms at a sam- 

pling rate of 40 MHz) 
HG high gain ( 'HG DOM' : one of the two DOMs in a tank) 

HLC 'Hard Local Coincidence' (in combination 'HLC hit' or 'HLC mode', describes an 

operation mode with full waveform readout) 
HV high voltage 

IC IceCube (in-ice part) 

IT IceTop 

LC local coincidence (of the two high-gain DOMs of a station) 

LG low gain ( 'LG DOM' : one of the two DOMs in a tank) 

MPE multiple photoelectrons (used in 'MPE discriminator', see also SPE) 

PE photoelectron (charge which a single photoelectron generates and is seen by the elec- 

tronics, including PMT gain, used as charge unit) 
PMT photomultiplier tube 

rms root-mean-square 

SLC 'soft local coincidence' (in combination 'SLC hit' or 'SLC mode', analogous to 

'HLC mode' to describe the operation mode; despite the name, 'SLC is not a 'coin- 
cidence') 

SPE single photoelectrons (used in 'SPE discriminator', see also MPE) 

UTC Coordinated Universal Time (reference time for all IceCube signals) 

VEM vertical equivalent muon (charge which a vertical muon generates in the IceTop DOM 

electronics of a tank, used as charge unit for IceTop signals) 



The primary mass determination from extended air showers (EAS) is notoriously difficult because the measure- 
ments are indirect and have to rely on models for the hadronization processes. Observables sensitive to the primary 
mass composition are mainly the height of the shower maximum (measured through fluorescence or Cherenkov 
emission) and the number of muons in a shower. The highest energy muons stemming from the first interactions 
in the higher atmosphere are most closely correlated to the mass of the primary nucleus. IceCube, in combination 
with IceTop, offers the unique possibility to observe these muons, typically with initial energies above 500 GeV, in 
the deep ice in coincidence with the mostly electromagnetically deposited shower energy measured at the surface. 
This provides a particularly powerful method for the determination of the mass composition. 

To mitigate the dependence on hadronization models, several alternative methods for studying mass compo- 
sition with IceTop are available. Other mass sensitive observables are for example: the shower absorption in the 
atmosphere at different zenith angles, the number of dominantly low-energy muons in the surface detector, and 
other shower properties such as shower age (state of shower development) and shower front curvature. 

The IceTop array has additionally been used to study high-px muons, PeV gamma rays and transient events, 
such as the radiation effects of solar flares. It also serves as a partial veto for the detection of downward-going 
neutrinos with IceCube and for direction calibration. 

This paper is organized as follows: First, we describe the technical design, deployment and commissioning 
of IceTop in Section [2] followed by the description of the front-end electronics in Section [3] and the trigger and 
data acquisition systems in Section [4] Next, the calibration procedures and the environmental conditions which 
influence the shower detection are discussed in Section [5] and [6] respectively. The next part deals with the data 
handling and analysis, beginning with the signal processing and data preparation in Section [7] and continuing with 
the description of the air shower reconstruction in Section [8] and the simulation of air shower events in Section [9] 
Finally the performance of the detector is demonstrated in Section 10 with examples from finished and ongoing 
analyses. The paper ends with a summary and outlook. In all descriptions we concentrate on the current status 
of the detector and the software (spring 2012) and refer to previous configurations only if there were significant 
changes which matter for analyses using those data. Acronyms used in the following text are explained in Table [T] 
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2. Technical design and deployment 



2.1. Design considerations 

The IceTop air shower array was proposed as an addition to the IceCube neutrino telescope to extend the 
telescope's capabilities for cosmic ray physics and to partially veto the background of down-going muons. In order 
to realize this addition in the most cost effective way, technical developments and infrastructure of IceCube were 
employed as much as possible. These considerations defined the following guidelines for the design of IceTop: 

• IceTop uses the same cables as laid out in trenches for the IceCube holes. That constrains IceTop's array 
size to 1 km 2 and the grid spacing to 125 m. Both together define the energy range from some 100 TeV to 
about 1 EeV. 

• Whenever possible IceTop uses the same detector hardware, electronic readout, triggering scheme and data 
acquisition as in the deep-ice part of IceCube (henceforth referred to as 'in-ice IceCube'). 

Using the same Digital Optical Modules (DOMs) as in the deep ice for detection of Cherenkov light with the same 
readout scheme leads to the choice of ice as detector medium. The ice is produced in tanks similar to the water- 
filled tanks used in the Haverah Park experiment JU and the Pierre Auger Observatory Q. Such ice tanks are well 
suited for the detection of shower particles because a particle produces enough light so that signal fluctuations are 
not dominated by photon statistics. For example, in the finally produced tanks a vertically through-going muon 
produces approximately 125 photoelectrons. The choice of ice Cherenkov detectors saved design and development 
costs and also simplified integrating the IceTop signals into the in-ice data acquisition system. 

Given the constraints on size, spacing, detector technology and basic components the following further consid- 
erations aim at optimization for air shower physics: 

• Air showers have a huge range of different energy depositions in the tanks. To increase the dynamic range 



offered by the DOM electronics (see Section 3.1 1 each tank has two light detectors running at different gains. 
Having two DOMs in each tank also allows for an experimental determination of internal signal fluctuations 
by setting both DOMs to the same gain during special calibration runs. 

Each detector station consists of two tanks for the following reasons: 

- By requiring events to trigger both tanks of a station ('local coincidence') signals from extensive air 
showers can be distinguished from the high (typically 2 kHz) individual event rate generated mostly by 
low energy showers. 

- Signals seen in both tanks of a station can be compared on an event-by-event basis to give a measure 
of the intrinsic physical fluctuations in the shower front, both in timing and in amplitude. 

- IceTop can be divided into two very similar sub-arrays, each comprising one of the two tanks of each 
station, to measure the fluctuations of reconstruction parameters, such as shower front curvature, lateral 
distribution, core location accuracy, angular resolution. 

- Single-station events that trigger both tanks of a single, interior station without triggering adjacent 
stations could be exploited to extend the IceTop reach to lower energies. Such events have primary 
energies around 30 TeV and have typically at most one muon with sufficient energy to reach the deep 
detector. 

The size of the tanks and their wall reflectivity were optimized for large light output and short signal pulses. 
Both, large tank size and high reflectivity, yield high signals but long pulse lengths (in the final tank design 
the decay constant is about 30 ns). 



2.2. General detector layout 

The IceTop air shower array is located above the in-ice IceCube detector at a height of 2835 m above sea level, 
corresponding to an atmospheric depth of about 680 g/cm 2 . It consists of 162 ice Cherenkov tanks, placed at 81 
stations and distributed over an area of 1 km 2 on a grid with mean spacing of 125 m. Figure [5] shows a plan view of 
the final IceTop array. In the center of the array, three stations (numbered 79, 80, 81) are installed at intermediate 
positions. Together with the neighboring stations (numbered 26, 27, 36, 37, 46) they form an in-fill array for denser 
shower sampling. Each station comprises two cylindrical tanks, A and B, 10 m apart from each other. The tanks 
are embedded into the snow so that their top surface is level with the surrounding snow to minimize temperature 
variations and accumulation of drifting snow. 

IceTop, as well as in-ice IceCube, took data already during the construction phase; the periods during which 
certain IceCube-IceTop configurations had been commissioned are reported in Table|2] 
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Table 2: List of the years when a certain configuration of IceCube (IC), IceTop (IT) became operational together with the numbers of the 
stations which were added in that year (fourth column). The numbers for IC include also the DeepCore strings. We will use abbreviations like 
IC79/IT73 for the constellation in 2010, for example. The fourth and following columns list stations using the numbers given in Fig. [2] The 
fifth column lists the tanks which have Tyvek liners (the others are coated with zirconia liner). However only the tanks installed in 2005 have 
higher reflectivities than the other tanks. The sixth column reports the DOMs which have the 'old' transformers with 43 CI and a short time 
constant, see Sections [3.2||3.5| and Table|4] The last two digits in this column (62 or 64) are the DOM numbers in a station according to Fig. [6] 



Year 


IC 

strings 


IT 

stations 


new IT stations 


Tyvek liner 


old transformer 
(43 Q) 


2005 


1 


4 


21 29 30 39 


all 


all 


2006 


9 


16 


38 40 47-50 57-59 66 67 74 




all 


2007 


22 


26 


46 55 56 64 65 71-73 77 78 




(77,46,56)-62 












(77,46,7 1)-64 


2008 


40 


40 


44 45 52-54 60-63 68-70 75 76 




(53,55)-62 












(52,55,62,68-70)-64 


2009 


59 


59 


2-6 9-13 17-20 26-28 36 37 




37-64 


2010 


79 


73 


8 15 16 23-25 32-35 41-43 51 




23-62 


2011 


86 


81 


1 7 14 22 31 79-81 


79 A, 80A, 81A/B 





600 



400 



200 



£ 



-200 



-400 
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-600 



-400 



-200 





X[m] 



200 



400 



600 



Figure 2: Locations of IceCube string holes and IceTop tanks with the IceCube Lab (ICL) in the center in the final configuration after 2010. 
The holes 81 to 86, belonging to DeepCore and not related to IceTop tanks, are not shown. IceTop stations are located next to IceCube strings 
(except for the in-fill station 81) and consist of two tanks, A and B. The irregularity of the array arises because tank locations were constrained 
by surface cabling and IceCube drilling operations. A denser in-fill array is formed by the stations 26, 27, 36, 37, 46, 79, 80, 81. 
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182 cm 



Figure 3: Tank dimensions in a cross-sectional view. Wall thicknesses are not to scale (see Table[5J. 
Table 3: Dimensions of IceTop tanks and assembly. The numbers given in units of m have tolerances of the order of 1 cm. 



component 


quantity 


value 


polyethylene tank 


height 


1.10m 




wooden extension 


0.20 m 




inner diameter 


1.82m 




wall thickness 


6 mm 




zirconia liner 


4 mm 


ice 


height 


0.90 m 


DOMs 


distance between centers 


0.58 m 


perlite 


thickness 


0.40 m 


outside tank 


polystyrene below tank 


(100 ±2) mm 




polyurethane foam around tank 


(50- 100) mm 



2.3. Tank design 

A schematic cross section of the IceTop detectors is shown in Fig.[3](dimensions are reported in Table|3]l. The 
tanks are made of black, cross-linked polyethylene, 6 mm thick, 1.1m high, with a 1.82 m inner diameter and are 
filled with ice to a height of 0.90 m. Most of the tanks have an integral diffusely reflective white liner made by 
dispersing zirconium dioxide powde^J referred to as zirconia, into High Density Polyethylene (HDPE) by extru- 
sioij^] The polyethylene, containing 6% zirconia by volume, is milled to fine white powde^jwith average particle 
size 45 fim. The tanks are produced by a rotational molding techniqu^] starting with molding the 6 mm thick, 
black outside layer and then covering it using the white powder to form a 4 mm thick layer. A final curing process 
leads to cross-linked bonds in the HDPE polymer structure which strengthens the tanks and finally determines the 
diffusive reflectivity of the liner. Eight tanks commissioned in 2005 and four tanks deployed in 201 1 have Tyvelj^] 
linings in form of bags loosely covering the tank walls (Table|2]l. 

Figure |4] shows measurements of the diffuse reflectivity of tank liners as a function of the wavelength. The 
measurements were done outside of the tanks and usually in dry condition. For Tyvek also a comparison of 
reflectivities of dry and wet material from the 2005 tanks is shown in the plot. The difference between dry and 



Supplier: Stanford Materials, Irvine, CA 92618 U.S.A. 
6 Manufacturer: PlastiScience, LLC, Smyrna, DE 19977 
'Manufacturer: Power King, Texas 
8 Manufacturer: PolyProcessing, Winchester VA 

'Dupont brand of a synthetic textile made of high-density polyethylene fibers 



8 



1 




Figure 4: Measurements of the wavelength dependence of diffuse reflectivity on samples from Tyvek bags ('uncoated' refers to the 2005 tanks) 
from tanks with zirconia liners. For more details see text. 



wet Tyvek is relatively large because the water fills the Tyvek matrix which is for dry material filled by air. The 
reflectivity depends on the ratio of the refractive indices of filler and matrix which is difficult to predict for Tyvek 
under real conditions in ice-filled tanks. The Tyvek tanks commissioned in 2011 (Table [2) exhibit in ice optical 
properties which are resembling much more the properties of zirconia tanks according to pulse decay time and 
signal size (see Table 15). The reason for this unexpected behavior could not be fully traced back because several 
other conditions were different for the first tanks installed in 2005. For example, in 2005 the water came from 
the station, while later the water was prepared in the drill camp and additionally filtered. For zirconia tanks all 
reflectivity curves are quite similar as demonstrated by the two examples in Fig. |4] In contrast to the Tyvek liner 
the zirconia liner does not absorb water. However, the real conditions for the liner-ice boundary are not known 
in both cases. In the simulation (described in Section [9.2. 1[ ) the measured reflectivity curves are scaled to match 
the observed optical properties of the tanks, which are the pulse height and decay time (for typical single-particle 
pulses as measured with muons, see Section 3.3 1. 

Figure [5] shows the tank assembly. Attached to the top of the tank is a wooden support structure (shown in 
brown in Fig. [5} that consists of 2" x 8" joists covered by hinged, plywood lids. The joists support the DOMs 
before the water freezes and the whole top serves as a structural platform to protect the detector after deployment. 
The hinged panels allow access for installing DOMs and let heat escape during the freezing process. The space 
between the ice surface and the lid is filled with expanded perlite (amorphous volcanic glass, expanded to low 
density with grain sizes of the order of 1mm [9|) for thermal insulation and light protection (see Table [9] for 
composition). 

The tank rests on a base of 10 cm thick rigid polystyrene insulation placed on top of a wooden pallet. The 
tank is mechanically strapped to the pallets with ropes. An enclosure for the Freeze Control Unit, described below, 
is mounted on the tank. On the right side of the Freeze Control Unit the drawing shows a tube into which poles 
carrying survey markers are inserted. Two vertical lines of temperature sensors are installed on the outer wall of 
the tank (not shown in Fig. |5) to monitor the freezing progress, one next to the Freeze Control Unit and one on 
the opposite side. There are 8 sensors per line in distances of about 15 cm, starting very near the tank bottom 
and ending about 15 cm above the water/ice level (measurements with these sensors are displayed in Fig. [9). The 
whole tank structure is integrated by applying expanding polyurethane foam around the tank from the pallet to the 
lid. The frame structure above the tank is only attached during the freezing process and carries sunshades to keep 
direct sunshine off the ice surface. 

The assembly took place in the US and the completed units were transported to the South Pole. The sunshade 
frames were installed before deployment and removed after freezing was complete. 



2.4. Deployment and commissioning 

IceTop deployment was integrated with the in-ice drilling and deployment. Tank locations were laid out during 
the survey for the surface cables of IceCube. Stations consisting of two tanks, 10 m apart, were placed nominally 
25 m from the top of the corresponding in-ice hole. Looking from the hole towards the tanks the left tank was 
labeled A' and the right one 'B' (Fig. [6), with the exception of station 48 where it was accidentally reversed. 
Tanks were set into trenches so the tops of the tanks were even with the snow surface to minimize drifting. After 
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Figure 5: Design rendering of an IceTop tank assembly before the outer foam insulation is applied. 

to IceCube Lab 




to IceCube Hole 



Figure 6: Sketch of the orientation of tanks and DOMs. Indicated is the position of the tanks relative to the nearest IceCube hole. Between both 
tanks is the 'surface junction box' which connects the 64 DOMs of the IceCube string and the IceTop tanks to the cable running to the IceCube 
Lab. The string DOMs have the numbers 1 to 60, the IceTop DMs from 61 to 64. The IceTop DOM numbering scheme is shown in the sketch. 
The DOMs 61 and 63 have high gains and 62 and 64 low gains (see also Fig.|15). 



the tanks were in position, the DOMs were mounted on hangers attached to the structural joists of the tanks (Fig.[5]l 
with the PMT half facing down. The equator of the DOM was set at 90 cm above the bottom of the tank so that the 
photocathode portion of the DOM was submerged. The two DOMs in each tank were placed along a diameter of 
the tank symmetrically about its center with a separation of 58 cm between the centers of the two DOMs. Figure|7] 
shows a tank with its two DOMs installed in the ice just before the tank was filled with perlite insulation and 
closed. 

A 'surface junction box' was placed between the two tanks of each station (Fig. [6]) and the four IceTop DOMs 
as well as the corresponding in-ice cable with wires for its 60 in-ice DOMs were connected to it. The connection to 
the IceCube Lab (ICL) at the center of the array was made by a surface cable carrying the wires for the 64 DOMs 
at each hole/station. The order of deployment of the IceTop stations followed the trenching and laying of surface 
cables. The surface cable also carries wires to provide power and housekeeping signals for the Freeze Control 
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Figure 7: Photograph of a frozen tank ready to be closed. The degasser unit, which remains frozen at the bottom, is visible under 90 cm of ice 
in the lower left. The fuzziness in the outline of the degasser is mainly the result of imperfections in the surface of the ice. 



Unit. At the end of the season, these wires were secured to the survey marker of each tank and are available for 
possible future use. 

Since the Freeze Control Units had to be kept warm, the tanks could not be filled until the surface cables were 
connected in the IceCube Lab. When the surface cables were laid and connected to power, the Freeze Control 
Units were installed and immediately powered up. Soon after, the tanks were filled with water from the in-ice hot 
water drilling system (except for the tanks commissioned in 2005, see remark in Section 2.3 I. At the beginning of 
the drilling season the water contained as much as 4000 ppm of propylene glycol, which remained from the system 
protection for the winter. A filter assembly reduced the glycol content to a level of less than 10 ppm in the IceTop 
tanks. The filtered water was transported to each station in an insulated water reservoir mounted on a dedicated 
sled. The water was filled to a depth of 90 cm and the surface was maintained at that level until the freezing was 
complete. 

Immediately after the tanks were filled, the Freeze Control Units were turned on with the tank lids closed 
for an initial two to three day period of cooling and degassing. During this period, sunshades were mounted 
on pre-installed frames that extended above the tanks (Fig. BJ. The tanks were then opened and the freeze-in 
process lasting about 50 days began. The surface of the ice in the tanks had to be kept clear of snow to prevent an 
accumulation that would have insulated the tank and slowed or stopped the freeze process. After an initial period 
of operating the degasser nearly continuously, the duty cycle was decreased in steps to the minimal level needed to 
insure clear ice. After the freeze was complete the top of each tank (40 cm between the ice surface and the lid) was 
filled with expanded perlite insulation and the top closed. The sunshades were removed, and the DOMs could be 
turned on and commissioned. The snow was backfilled and groomed around the tanks. The Freeze Control Units 
were retrieved from the tanks and stored for future use. 



2.5. Water freezing 

A natural process that allows the water to freeze from the top to the bottom was adapted to make clear, crack- 
free ice. When ice freezes in a confined volume from the top down, two problems must be managed. As the ice 
front advances, air is trapped in the remaining water and would soon rise to saturation level and make bubbles if it 
were not removed. The expansion at the phase transition from water to ice also has to be dealt with. 

The principle of the freezing operation is sketched in Fig. [8] The air is removed using a commercially available 
contactor (A in Fig. [8} consisting of a cluster of tiny tubes made out of a semi-permeable membrane passing 
through a plastic cylinder. A submersible pond pump (B in Fig. [8} with a five-micron filter to keep any sediment 
from fouling, forces water through these tubes. The contactor is kept at a modest vacuum (130hPa). About half 
the dissolved air is removed in one pass with this setup. Except for the vacuum line and the power supplying 
the pump, the unit is self-contained. A negative pressure ensures any leaks are inward, and the resulting water is 
removed by a water tolerant vacuum pump. Tests showed that it is sufficient to run this assembly at the bottom of 
the tank. If the tank were perfectly insulated, the densest water (4 °C) would concentrate at the bottom and only 
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(a) Gas exchange contactor 
® Water circulation pump 



Figure 8: Schematic of the operation principle of the Freeze Control Unit: Initially the tank is filled until the water level spills over to the 
second container with the sump pump which pumps overflow water to the outside (third pipe). At the bottom of the tank, a pump pushes water 
through the degasser unit, keeping the water circulating until all the water is frozen. 



the lower portion of the tank would be degassed. Even though the sides of the tank are insulated, there is enough 
thermal transfer so that a stable temperature gradient from top to bottom was never set up. The small amount of 
heat loss through the sides is sufficient to keep the water in the tank well mixed during the freezing process. The 
temperature of the entire volume, which remains near 0°C, and the dissolved gas content are nearly independent 
of depth. With the pump operating at full capacity of 2 liter per minute, the entire 2500 liter volume of water in the 
tank can be processed about once per day. 

Removal of the water as the ice front advances is also straightforward. The ice rapidly bonds to the sides 
of the tank with sufficient strength to form a plug. The pressure builds only high enough to force the water out 
through a pipe (Fig. [8]) open at the bottom of the tank that extends into the insulated and heated Freeze Control 
Unit enclosure. The short horizontal section sets the initial water level in the tank. Inside the Freeze Control Unit, 
the water collects in a sump which is pumped out when full. Pumping water out rapidly when needed prevents the 
water from freezing and plugging that discharge hose which is exposed to about -30 °C temperatures. The Freeze 
Control Unit controls the circulation and vacuum pumps for the degasser, and the sump pump. Resistive heaters 
are used to keep the temperature inside the Freeze Control Unit at 30 °C. 

Freeze history. The history of a freezing tank is illustrated in Fig. [9] depicting the depth profile of water/ice tem- 
peratures in the tank during the freezing process. The temperatures are measured through the line of temperature 
sensors mounted vertically on the outside. Initially the tank is empty so all of the sensors measure the ambient air 
temperature (here around -20 °C). Water is then added to the tank and all but one of the sensors now follow the 
water temperature. The exception is one sensor that is located above the surface level of the water. At this stage the 
tank lid is kept closed so the water can cool uniformly and the degassing unit can reduce the concentration of air 
in the water (see description above). When the water has reached equilibrium and is near the freezing point the lid 
is opened, causing an immediate drop in the 'free air' sensor (the highest in the sensor line) to near ambient level. 
The water begins to freeze from the top down. As the freezing front passes each sensor in turn the temperature 
drops below °C. The thermal profile of the sensors near the surface follows the ambient temperature more closely 
than that of the deeper sensors as would be expected from the need for the heat to diffuse through the block of 
ice. Eventually all of the sensors fall below zero, but the thermal gradient remains high since there still is liquid 
water in the center of the tank even though the sides are completely frozen. At this stage water is also still being 
drained into the sump, further indication that the freeze is not complete. Finally, when the last water freezes the 
temperature gradient collapses and water evolution ceases. The freeze control unit is now disconnected, the perlite 
insulation placed on top of the tank, and the lid screwed shut. At this point the 'free air' sensor now joins the others 
as it is buried in the insulation and thermally better connected to the ice than it is to the outside air. 
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Figure 9: History of the freezing process as recorded by the temperature sensors which are mounted to the wall of the tank at different heights 
(8 sensors starting from the tank bottom to 15 cm above the water/ice level, see Section [23) . In the plot the lowest curve, staying at about 
-20 °C, is measured by the highest sensor (above the water), the next curves following upwards are measured by the sequentially lower sensors. 



2.6. Survey of the tanks 

Surveyors determined the location of the center of each IceTop tank and the location of the survey marker on 
the tank. From those two measurements and the known tank parameters, the location of the center of each IceTop 
DOM is determined to an accuracy of 5 cm. The tops of the deep holes of IceCube were determined to a similar 
accuracy in the same coordinate system. Initial location of the deep DOMs was determined from the cable payout 
records and pressure sensor data. In a second stage LED light pulses generated in the DOMs (flashers) were used 
to determine relative vertical offsets of strings. Reconstructed muon data and flasher data are used to monitor for 
shearing or other changes over time, which have so far not been observed. Locations of deep DOMs are known to 
an accuracy of better than 1 m. 

The IceTop tanks lie roughly on a tilted plane with the top of the tanks at z — 1944.54 m for the lowest point 
(tank 18A in Fig. [2} and at z = 1950.08 m for the highest (tank 74A) where z — is defined in the center of 
IceCube. 



3. Front-end electronics 

The pivotal components of the detector are the Digital Optical Modules (DOMs). They detect the Cherenkov 
light produced by the shower particles in the ice of the tanks and convert these signals into digital information, 
which is sent to the central counting house, the IceCube Lab. 

IceTop DOMs have the same hardware as all other DOMs in IceCube, but some details of the data acquisition 
system are different due the different physics requirements and different environmental conditions. For example, 
the pulses generated by air showers deliver in general more charge and thus the pulse thresholds are higher; another 
difference is that the IceTop DOMs are subject to the seasonal temperature changes while the temperature in the 
deep ice is extremely stable. 



3.1. Digital Optical Module (DOM) 

Each DOM 1 10 1 consists of a 10" Hamamatsu photomultiplier tube (PMT) IfTTII and electronic circuitry, con- 



tained inside a 33 cm glass pressure sphere, as shown in Fig. 10 a. The PMT is shielded from Earth's magnetic 
field by a mu-metal grid. The two DOMs of a tank are operated with different PMT gains, high gain (HG) and low 
gain (LG), to adapt the dynamical range to the extremely different signals in air showers (TableQ. 

Data taking, triggering, digitization and communication with the IceCube Lab are controlled by an FPGA 
(Field Programmable Gate Array) with an on-chip CPU on the DOM mainboard. The CPU permits remote DOM 
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Figure 10: The IceCube Digital Optical Module (DOM): a) Schematic view of a DOM. b) Analog signal part of the block diagram of the DOM 
mainboard. 



firmware updates. Timing of the DOM is controlled by a free-running 20 MHz oscillator which is regularly cali- 
brated with the master clock in the IceCube Lab. In addition, each DOM mainboard is equipped with an LED for 



calibration, see Section 5.1 (there are other LEDs in the DOM which are used to study the optical properties of the 
deep ice but are not used in IceTop). 

A detailed description of the Digital Optical Modules can be found in [ 10 1 . Some characteristic parameters for 
the IceTop DOM operation are summarized in Table [4] 

3.2. Signal capture and digitization 

The PMT signals are captured by the DOM electronics when the peak pulse voltage surpasses a programmable 
discriminator threshold. The relation between the pulse voltage, which the discriminator triggers on, and the pulse 
charge, which is the calibration quantity depends on the pulse shape and thus on the front-end electronics. Because 
of the importance for IceTop signal and threshold calibration we briefly describe the two electronic paths for pulse 
discrimination and recording. 

The PMT's cathode is kept grounded with the high voltage (HV) applied to the anode from where the signals 
are read out. The signals are decoupled from the HV by a wide-band inductive transformer (also referred to as 
'toroid') which is expected to have a higher reliability as compared to the more common capacitive decoupling. 



The effect of the droop caused by a transformer and its correction will be discussed below (Section 3.5 1 



The signal paths described here can be followed in the schematic given in Fig. 10 b. From the PMT's anode the 
signal passes the transformer and is then split into a path to the discriminator trigger and two paths to the signal 



recording electronics. A 'Pulser', feeding into the same paths, is used in the calibration procedure (Section 5.1 1. 

If the signal passes the discriminator threshold, the PMT output is sampled by a custom-made integrated 
circuit called 'Analog Transient Waveform Digitizer' (ATWD) lfT3l in 3.33 ns wide bins for a total of 128 bins, 
corresponding to a total sampling time of about 427 ns. The signals on the ATWD path are delayed by 75 ns to 
allow for the trigger decision. Each ATWD has four channels, labeled ATWD0 to ATWD3, three of which are 
used for waveform recording. The fourth channel is not used for regular data taking but, depending on the running 



mode, other information can be recorded using a multiplexer (MUX in Fig. 10 b), see details in IfTUl . To minimize 
dead-time, each DOM is equipped with two ATWD chips which are used alternately. 

The inputs to the channels 0, 1 and 2 are amplified by nominal gain factors 16, 2 and 0.25, respectively (Table 
|5j. In the case of a discriminator trigger, the FPGA opens the input channels of one of the two ATWD chips and 
the three scaled signals are recorded synchronously. After the analog sampling, 128 Wilkinson ADCs digitize with 
10 bit precision each of the 128 ATWD cells of a channel in parallel. The three channels are digitized sequentially, 
starting with the highest gain channel; channels with lower gain are only digitized if in the higher gain channels 
any bin is above 768 counts. This takes about 30 jus per channel and up to 100 /is for all three channels. The 
three different gain channels of the DOMs make up for an effective dynamic range of 16 bit. Measuring charges 
in units of PE, the charge yield of a single photoelectron, the effective dynamic range of the high-gain DOM in a 
'standard' tank (Table [6]) spans from about 0.15 PE to about 10000 PE. The conversion from mV to PE is done 
by using the numbers in Table [6] (in this case for the 'standard' tank) which are specific to the multi-PE pulses 
seen in IceTop tanks and not to the pulse shape from single photoelectrons. The standard high-gain IceTop pulse 
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Table 4: Important parameters of front-end components. 



component quantity 



value 



PMT 



transformer 



quantum efficiency at 390 nm (near maximum) 
quantum efficiency, av. for 300 - 650 nm 
cathode area 
high voltage, gain 

saturation current (50% from linearity) 

time constant (at -30 °C): 
front-end impedance 



~ 25% HH 

~ 18% HQ 

~ 550 cm 2 |[T2l 

HG:5xl0 6 at~1250V 

LG: 10 5 at -750 V 

HG: 140 mA 

LG: 50 mA 

1 — 2/x s (old), 12 yus (new) 
50 Q (old), 43 Q (new) 



Table 5: Characteristics of the three ATWD readout channels (nominal or approximate values). The smallest digitization value of each channel 
is given in the third column as LSB = least significant bit. 

channel gain LSB [mV] saturation [mV] 

(high gain) ~ 16 0.125 -100 

1 (medium gain) 2 1.0 -800 

2 (low gain) 0.25 8.0 -7500 



Table 6: IceTop pulse characteristics. The numbers are average or approximate values, actual values can deviate. The notations 'old' and 
'standard' tanks refer to the tanks commissioned in 2005 and after 2005, respectively. 



tank liner 


'standard' tanks 


'old' tanks 


PMT gain 


HG 


LG 


HG 


LG 


rise time (10-90%) 


10 


10 


10 


10 


decay time (observed) [ns] 


31.0 


31.0 


43.5 


43.5 


decay time (input for simulation) [ns] 


26.5 


26.5 


41.9 


41.9 


VEM charge [PE] 


125 


110 


200 


180 


fraction of charge in peak bin 


0.07 


0.07 


0.05 


0.05 


peak voltage of 1 VEM [mV] 


110 


2.2 


100 


2.0 


peak voltage of 1 PE [mV] 


0.88 


0.02 


0.5 


0.01 


total charge of 1 PE [mV] (sum of ATWD bins) 


12.6 


0.29 


10.0 


0.2 


saturation in PE 


8000 


125 000 


14000 


225 000 


saturation in VEM 


64 


1140 


70 


1250 



has a peak voltage of 0.88 mV/PE corresponding to about 7% of the total charge in the peak bin (for the single 
photoelectron pulse shape these numbers are about 40% or 5 mV/PE in the peak bin). 

In response to a trigger, the ATWDs are launched in discrete time steps of 25 ns, phase synchronized to the 
frequency-doubled 20 MHz DOM clock. The final time resolution is obtained from the position of the waveform 
in the ATWD with respect to the start time. The pulse time determination is described in Section|7] 

In parallel, the pulses are digitized continuously by a commercial 10-bit fast ADC (fADC in Fig. 10 b) at a 
sampling rate of 40 MHz. The fADC inputs are shaped with a time constant of 180 ns limiting the bandwidth to 
adjust to the lower sampling speed. In response to a trigger (Section|4]i the downstream DOM electronics records 
fADC data for 6.4^ s. The fADC data have not yet been used in IceTop analyses but could be used for studies of 
signal tails extended in time. 

If DOM local coincidence or certain other conditions as descibed in Section |4] are satisfied, the ATWD and 
fADC digitization streams are losslessly compressed and written to memory. If a discriminator has triggered and 
an ATWD is available, at least the sum of the ATWD bin contents, called chargestamp, and a timestamp is extracted 
and written to memory. The ATWD recording is not possible while a waveform acquisition is still in progress or 
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Figure 11: Typical IceTop waveforms. The most frequent shape has a steep rise, a single peak and an exponential tail. The waveform is here 
not droop corrected, but in panel a) the undershoot below the baseline after the pulse is very small, a feature of the 'new' transformers installed 
in DOM 09-61 (Table[3j. An example for the 'old' transformer characteristics is shown in panel b) with a quite distinct undershoot. This DOM 
21-63 belongs to the tanks installed in 2005 (Table[3} which exhibit longer decay times of the waveforms (Table|6j. In the left plot the blue line 
depicts the slope of the leading edge with the extrapolation to the baseline defining the pulse time (red circle), see Section|7] 



when both ATWDs are busy (see also the discussion on dead time in Section 4.2.4 1. 



3.3. Pulse shapes and discriminator thresholds 

As described above, the readout of signals is initiated if a pulse crosses a voltage threshold set by the DOM 



discriminator (upper path in the schematic in Fig. 10 b). Any relation between voltages and charges depends on the 
pulse shape, which is determined by the tank properties and the electronics. Therefore we first discuss the typical 
shapes of tank pulses before turning to discriminator thresholds. 

3.3.1. Pulse shapes and VEM units 

Most of the pulses have a similar shape given mainly by the reflectivity of the tank walls because the spread 



of crossing times of the shower particles are short compared to the decay time due to reflections (Fig. Hi. They 
can mostly be described by a superposition of standard pulses which is well represented by the signal of a single 
muon. The signal charge of a single muon is also used as a calibration standard for normalizing the signals of 
different tanks. The normalized charges are expressed in units of 'vertical equivalent muons' (VEM) determined 
by calibrating each DOM with muons (see Section[5]). 

The standard pulses are characterized by a fast rise of the leading edge, about 10 ns from 10% to 90%, and an 
exponential decay with a time constant of on average 31 ns for most of the tanks (Fig. 11a) and 43.5 ns (Fig.[TTjb) 



for the 8 tanks installed in 2005 (Section [2. 3| >. If not otherwise mentioned we use as nominal pulse shapes those 
obtained for the tanks which were installed after 2005 (Table|6]). In Table|6]two decay times are given: as observed 
and as put into the simulation. The difference is due to the electronic shaping which is applied also in the simulation 
as described in the next section. 

3.3.2. Pulse shaping by the electronics 

IceTop analyses are particularly dependent on a good simulation of pulse shapes to correctly reproduce the 
experimental thresholds: Unlike the case for in-ice DOMs where the thresholds are set on fractions of a PE, the 
multi-PE thresholds in IceTop are more affected by the actual pulse shape generated by many PE's. 

For the simulation, pulse shaping is described by transfer functions for the different paths with which the PMT 
output has to be convoluted. The software simulating the DOMs [ 14 1 currently uses the following transfer function 
for the three ATWD channels and the discriminator path: 

f(t) = —=exp\ ^—n (1) 



t V^TT \ 2<X 2 

with the parameters <x and t m characterizing the general shape of the signal. The parameter <x depends on the 
transformer version (0.41 for the 'old', 0.35 for the 'new' transformers) and t m is set for both to 13.53 ns. The 
function has a maximum at the time t„ mx which is related to the parameter t m by: 

tmax tm & ■ (2) 
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Figure 12: Comparison of the transfer function currently used in the detector simulation, as described by the function JTJ with the parameters 
given in the text, with the measured transfer functions of the ATWD and discriminator paths (for 'new' transformers). 




Figure 13: Threshold monitoring plots for the high-gain DOM in tank 26A. a) The peak voltage of the first arriving pulse of the waveform 
shows a sharp rise at a threshold of 24 mV. b) The threshold in terms of charge is smeared out. The charge threshold can be defined by 50% 
point of the leading slope of the charge distribution as shown on the plot. 



A detailed study [ 15] found that the ATWD and discriminator paths have similar, but not identical, transfer func- 
tions. Both can be described by |Q}, but the parameters <x and t m are significantly different resulting in narrower 
pulses at the discriminator than at the ATWD. The shapes of the transfer functions with the measured parame- 



ters for both paths and those currently used in the simulation are displayed in Fig. 12 After convolution of the 
PMT pulses with the transfer functions, differences between both paths are reduced and a single set of simulation 
parameters gives a reasonable approximation. 



3.3.3. Discriminator thresholds 

Each DOM has two voltage sensitive discriminators, called SPE (single-PE) and MPE (multiple-PE), which are 
generally used for triggering waveform capture. Both discriminators feature programmable thresholds with ranges 
differing by a factor 10. In IceTop, both discriminators are used for multiple-PE thresholds. In high-gain DOMs, the 
MPE discriminators are used for triggering on air showers with thresholds of about 20 mV corresponding to a signal 
charge of about 23 PE for a standard pulse (see Fig. [T3j ), while the low-gain DOMs use the SPE discriminators for 
the same purpose with thresholds of about 4mV corresponding to 270 PE |fl6l . With these thresholds, the single 
DOM rates are about 1600 Hz. 

The SPE discriminators of the high-gain DOMs are used to record 'scaler rates' at various thresholds (from 
about 0.5 PE to 30 PE) without triggering waveform capture. An example for the threshold dependence of rates is 
shown in Fig. [14] The scaler rates are monitored continuously for solar flares, gamma ray bursts and other transient 
events. 



3.4. PMT gains and saturation 

The two DOMs in each tank are operated at different PMT gains (TableQ: The high-gain DOMs are nominally 
set to a gain of 5 x 10 6 , half the standard gain for in-ice DOMs. At this gain the saturation current (deviation by 
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Figure 14: Count rate as a function of discriminator setting of different DOMs. Rates are averaged over the interval from 17 February 22:00 UT 
to 18 February 02:00, which defines the "base count rate" for a study of a Forbush decrease observed in February, 201 1 1 17 |. The line shows 
the expected values based on a simulation 1181 . 



50% from linearity) is about 140 mA corresponding to about 8 000 PE for the standard IceTop pulses (Table [6]). 
The low-gain DOMs are set to a gain of 10 5 and have a saturation current of 50 mA which corresponds to about 
125 000 PE for standard IceTop pulses. Although the gain of the low-gain DOMs is 50 times smaller than of the 
high-gain DOMs, the number of photoelectrons at saturation level increases only by a factor of 15. The available 
16 bit dynamic range of the ATWDs is about 10 000 PE for high-gain DOMs and 450 000 PE for low-gain DOMs 
(standard pulses), thus covering the non-saturated PMT range. The change from high-gain DOM to low-gain DOM 



signals is fixed to a much lower PE value, corresponding to less deviation from linearity (Section 5.2.5 1. Since the 
pulse shape stays very similar in a wide range (except for very low pulses and beyond saturation) the saturation 
can also be expressed in terms of VEM (for the 'standard' pulse, see Table |6j. The low-gain DOM saturation 
determines the IceTop tank saturation which is typically of the order of 1000 VEM. 



The implementation of saturation effects in the simulation is described in Section 9.2 



3.5. Droop 

The transformer described in Section 3.2 leads to a characteristic high-pass frequency response for all signals. 
While the time constant is large compared to IceTop pulse widths, the later parts of each waveform show a 'droop' 
in amplitude, followed by undershoot below the baseline after the pulse ends. The pedestals of the ATWD and 
fADC are set to about 10% of the maximum scale, to permit waveforms with below-baseline excursions (that 
effectively also reduces the available 10 bit range). 

Choosing the time constants of the transformer requires a compromise between loss of precision in high- 
frequency signals, and large droop after large signal amplitudes. The first DOMs produced had transformers with 
short time constants of about 1 fis at -30 °C leading to a strong droop effect (Fig. [TT|b), whereas DOMs produced 



later (since 2005) have larger time constants of about 12/vs at -30 °C with reduced droop (Fig. 11 a). In IceTop, 
79 of the 324 DOMs or 27% of all DOMs have transformers with large droop, see Table|2] 

The droop effect can be described empirically by an impulse response with two time constants T12 Ifl9l , 

S( t ) S(t) - N ((1 - f) e- f/T ' + / e-'^ ) , (3) 

where is a normalization constant. The formula is used in a discretized form for simulation and corrections of 
measured pulses (see Section [TJ. The droop constants /, T\, and have been measured in the laboratory as a 
function of DOM temperature. The temperature dependence is (T in °C) 

B 

ti(T) = A + 



1 + e- T ' c (4) 
t 2 (T) = 0.75 n(T) 

Using a single-tau approximation by setting / = the parameters are on average: 

A = 0.4 yus, B = 5.0yus, C = 16 °C (old transformers) 
A = 1.0yus, fi = 50.0yus, C = 25°C (new transformers) 
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In the IceTop tanks the DOM mainboard temperatures vary between about -20 °C and -40 °C throughout the 
year (in the deep ice the temperatures are stable), yielding for the new transformers t values between about 9 
and 17 /is. This has a major impact on the signal extraction and requires a continuous monitoring of the DOM 
temperatures. The inversion of the formula ([3]) is used to remove the droop in the data, see Section|7] 



4. IceTop triggers and data acquisition 

4.1. General scheme 

In normal operation, the DOMs act as autonomous data collectors, initiating the capture and digitization of 
waveforms in response to one of many triggering conditions (generated either internally to DOM or from within 
a full IceTop station). A DOM is capable of sending a variety of data records, depending on selected operation 
modes, to the IceCube Lab for consideration in detector-wide triggers and event construction or for calibration 
and monitoring purposes. IceTop data included in these triggered detector- wide readouts, known as an 'event', are 
available to a suite of filter algorithms, selecting the most interesting events for immediate transmission north via 
a dedicated satellite link. 

Trigger scheme. To select physically interesting events and to save bandwidth for data transfer IceCube has a 
multi-level trigger system including both hardware and software decisions. In the following an overview of the 
different 'triggers' is given, which also serves to define the terminology: 

Discriminator Trigger: Normally the capture of the waveform by an ATWD is initiated by the FPGA, when a 



pulse passes a discriminator threshold; autonomous within a DOM; described in Section 3.3 



DOM Launch: The initiation of the process of capture and digitization of waveforms and data preparation for 
transmission organized by the FPGA according to the operation mode. Normally, a DOM is launched by the 
discriminator trigger, but the trigger can also be generated internally by the FPGA, for example for a baseline 
determination (Section |4.2.3| ). To be launched the DOM has to be in the 'ready' state (see the discussion of 
dead times below). 

Hit: The basic information recorded upon a DOM launch is referred to as a 'hit' . 

Local Coincidence (LC): Each DOM in IceTop is connected by a network of cable connections to neighboring 
DOMs in the same station. The presence of a logical level set by the discriminator trigger in neighboring 
DOMs, the so-called Local Coincidence condition, can be used to set the content of hits sent to the IceCube 
Lab. Described in Section |4.2.1| 

Software trigger: A series of software algorithms formed in the IceCube Lab using DOM Launches satisfying 
the LC condition determine when data from the entire IceCube detector should be recorded into events. 

Online filter: A subset of events are selected for immediate transfer north by satellite link using a set of filtering 
algorithms that select the physically most interesting events for later detailed analysis. All events are archived 
on magnetic tape storage. 

Operation modes. How much information a hit contains when being transferred to the IceCube Lab, is determined 
by the operation modes for the readout of DOMs. Two different modes were foreseen for IceCube [ 10 1: 

SLC mode: In the 'soft local coincidence' mode data from all launched DOMs are transferred to the IceCube 
Lab. The data format of a hit contains a header with coarse information about the hit and in addition the full 
waveform data for 'HLC hits', that are hits with a fulfilled LC condition ('LC tag'). For hits which have 
no LC tag ('SLC hits') only the header with the coarse information is transmitted. For IceTop the coarse 



information is a timestamp and the integrated ATWD charge. See Section 4.2.2 
HLC mode: In the 'hard local coincidence' mode only those DOMs transfer data which have an LC tag. 
In both modes additional conditions can lead to full waveform transfer, such as in case of the calibration triggers 



described below in Section 4.2.3 The baseline-operation of IceCube is the SLC mode 1 10] which, however, was 
not applied before 2009. For IceTop the standard SLC mode has been modified so that the ATWD digitization 
is not aborted if there is no LC tag because the integrated ATWD charge is always transmitted in the SLC mode 
(for in-ice DOMs only fADC information is used for the SLC transmission so that the ATWD digitization can be 
aborted if there is no LC tag). 
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Figure 1 5 : Schematic of the Local Coincidence configuration of high-gain (HG) and low-gain (LG) DOMs in an IceTop station. The coincidence 
requirements are described in the text. The DOM numbering is explained in Fig. [6] 



4.2. Decisions at detector level 
4.2.1. Local Coincidence 

For standard data taking the full-waveform readout of a DOM is initiated if a local coincidence (LC) of a 
station is fulfilled. This results in a single station trigger rate of about 30 Hz compared to about 1600 Hz of a single 
high-gain DOM at a threshold of about 0. 1 VEM. The single DOM rate is dominated by small showers with a good 
fraction caused by single muons passing a tank (see Section 5.2.2 1. 

To establish a local coincidence the high-gain and low-gain DOMs of a station are interconnected via cables 
as shown in Fig. 15 The following simplified description of the local coincidence decision focuses on the IceTop 
applications (details can be found in iflOll ). If a high-gain DOM passes the discriminator threshold a logical signal 
is sent to both DOMs in the other tank, which is used to define a 1 fis coincidence time window. In the other tank, 
a DOM of either gain is read out if it had a discriminator trigger within this time window. Thus, in the normal 
configuration, station data always consist of signals from both high-gain DOMs and optionally additional signals 
from low-gain DOMs. However, due to the fact that the LC condition is only based on the discriminator trigger, 
dead time (see below) can cause data from a high-gain DOM to be missing in the rare cases where both ATWDs 
are busy. 

The connections between high-gain and low-gain DOMs and an additional, normally unused connection be- 
tween the two low-gain DOMs (not shown in Fig. 15 i add also redundancy for the case that one or both of the 
high-gain DOMs stops functioning. In such an event low-gain DOMs can be switched to high-gain and an LC 
connection can be reestablished. 

A fulfilled LC condition will result in storing the digitized waveforms into the local DOM memory from where 
it will be transfered to the IceCube Lab as described in Section |4.3.1| 



4.2.2. Soft Local Coincidence 

As explained above, 'SLC mode' refers to an operation mode which reads out every DOM for which a signal 
passed the discriminator threshold, at minimum recording the sparse 'SLC hit' information of integrated charge and 
timestamp. Normally, only HLC hits include the full waveform information. Including SLC hits in the analysis 
is useful for detecting single muons in showers where the electromagnetic component has been absorbed (low 
energies, outer region of showers, inclined showers) as well as for vetoing air shower background for in-ice studies. 



4.2.3. Calibration launches 

Under certain conditions full waveform readout can be initiated without the LC condition being fulfilled. Cur- 
rently there are two such special DOM launches, both serve calibration or monitoring purposes. VEMCal launch: 
This launch of single high-gain DOMs collects single muon hits for calibration of the VEM unit. It is described 
in Section 5.2.2 Beacon launch: In order to determine the ATWD and fADC baselines from 'empty events' all 
DOMs are launched at a rate of 1 Hz without requiring that a discriminator trigger has fired. 



4.2.4. Dead time and unusual hit patterns 

In order to achieve a DOM launch the DOM must be in the "ready" state. Several conditions may cause a 
DOM to be not ready. If the DOM got a valid LC the DOM is in a "busy" state for 6.4 ^us (the time the fADC 
records the waveform). If the DOM did not get a valid LC a new DOM launch is disabled until the DOM knows 
it did not receive a valid LC. This time is currently set to 2.5 yus for IceTop. For a launch rate of 1600 Hz, this 
creates a dead time probability of 0.4%. After that time the DOM is reenabled if there is a free ATWD, but if both 
ATWDs are busy the DOM remains disabled until an ATWD becomes available. For IceTop the SLC mode requires 
always the ATWD digitization to be completed (see the definiton of the SLC mode above), which takes about 30 fis 
(Section 3.2 1. Assuming the 1600 Hz of launches is uncorrelated, there is about 5% probability of overlapping the 
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Figure 16: Schematic of the data flow and the software components of IceCube data acquisition system. The data from the DOMs are collected 
by the StringHub modules (residing in the DOMHub computers) which send data to the trigger modules. Upon a global trigger the Event 
Builder gathers the corresponding data of the trigger window and sends the data to the online processing and filtering system PnE 



processing from two DOM launches, and a dead time probability of about 2.5% that a discriminator crossing will 
not result in a DOM launch. When operating in HLC mode, if there is no LC condition, the ATWD launch can 
be aborted, in which case the ATWD in question can be re-enabled after about 6.5 ps (the LC decision time plus 
about 4 yt/s), and the LC-coincidence dead time is reduced to 0.5%. In HLC mode, other cases occur under unusual 
circumstances; for example, if the high-gain DOM in one tank (say, A) and the low-gain DOM in the other tank 
(B) are the only DOMs with a discriminator trigger, then a full report for the low-gain DOM will be generated, but 
no report will be made for the DOM launch in tank A. This might occur if a single low energy electron occurs in 
proximity to the low-gain DOM in tank B, so that the light signal in the high-gain DOM is below threshold. 

4.3. Data acquisition and processing 



An overview of the data acquisition system of all IceCube components is given in Fig. 16 In the following we 
explain the features relevant to IceTop operation. 

4.3.1. Data transfer to the IceCube Lab 



Each IceTop DOM is connected to a computer, called DOMHub (Fig. 16 1, by a single twisted pair of wires. 
This differs from the setup for in-ice DOMs, where the lower counting rate allows each twisted pair to carry the 
signals from two adjacent DOMs. The DOMHubs are standard personal computers equipped with PCI digital 
optical module readout (DOR) cards. Each DOMHub accommodates eight IceTop stations (32 DOMs). 

A fulfilled DOM launch condition, usually the LC condition, will result in storing the digitized waveforms of 
the highest-gain ATWD channel which is not saturated (the lowest-gain channel if all are saturated) and of the 
fADC into the local DOM memory. In regular intervals the DOMHubs request these data to be transferred. In SLC 
mode, time and integrated charge information is always transmitted and for HLC hits (and some special DOM 
launches, see Section |4.2.3| l in addition the full waveform is transmitted. 

4.3.2. Data processing in the IceCube Lab 

In the DOMHubs the data are prepared for further processing by the StringHub modules. First the timestamps 



of each hit, which are calibrated to a common time reference by the RAPcal system (Section 5.1 1, are converted 
to UTC time. The hits are then time sorted and grouped together within time windows of a few microseconds. 
Subsequently the HLC hits are sent to the trigger units of the in-ice detectors and of IceTop where trigger algorithms 
are applied (see below). The triggers define readout windows, which differ for the different sub-detectors and 
triggers. In IceTop the window for the normal air shower events is +10/^s before the first and after the last hits 
which produced the trigger. The IceTop calibration trigger has a window of ±1 ps. The global trigger merges 
all overlapping individual trigger windows for all sub-detectors. Based on this global readout window, the event 
builder collects all launches from all sub-detectors that happen during the global trigger window. 

In the common IceCube event builder the data selected by the global trigger are transformed into the DAQ data 
format, sent to the online processing and filtering (PnF) and are stored on tapes. The tapes can only be shipped to 
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Table 7: Summary of IceTop rates. 



name condition rate [Hz] 
hit rates 

DOM (single) DOM launch 1600 

HLC (all DOMs) LC fulfilled 3000 

triggers 

IceTopSMT > 6 HLC hits 30 

IceTopMinBias HLC prescale 10 4 0.3 

IceTopCalibration VEMCal launch 30 



the north during the Antarctic summer. Therefore part of the data, as defined by the filters, are sent via satellite 



transmission. In addition, the data for monitoring transient events (Section 4.3.5 i are also sent via satellite 



4.3.3. IceTop trigger 

The current basic trigger for IceTop air shower physics is the IceTop Simple Multiplicity Trigger (IceTopSMT) 
which requires 6 HLC hits (usually from 3 stations) in IceTop DOMs within 6 /is: 

IceTopSMT : > 6 HLC hits. (6) 

The readout window starts 10/vs before the trigger window and lasts until 10/vs after the last of the 6 hits. Upon 
an IceTop trigger the whole detector is read out and similarly, when an in-ice HLC trigger is formed, the previous 
interval is recorded to capture any activity in IceTop. 

In addition a 'minimum bias' trigger, IceTopMinBias, triggers on any HLC launch in IceTop with a 'prescale' 
factor of 10 4 , that means only every 10 4 th HLC launch triggers. 

All IceTop VEMCal launches (Sectio n [4.2. 3[ ) are selected by the trigger IceTopCalibration. The treatment of 
calibration triggers is described in Sectionp] 

A summary of various trigger rates is given in Table [7] 

4.3.4. IceTop online filter 

For satellite transmission the data have to be filtered for data reduction. The filters, summarized for IceTop in 
Table|8] select events according to various physics cases or for monitoring purposes. For some filter classes only a 
certain fraction of events is transmitted (prescaling). 

The basic IceTop filter. The most important IceTop filter, called IceTopSTA3, is based on events triggered by 
the simple multiplicity trigger, IceTopSMT, and requires that at least 3 stations have an HLC hit. During the 
construction phase events with less than 8 stations were prescaled for satellite transmission. From 2012 on all 
events with 3 or more stations are transmitted without prescale, though with partly only condensed information 
employing the so-called SuperDST format (see below). The three-station condition has an energy threshold of 
about 300 TeV. 

To allow adjusting to different transfer and/or analysis modes the events passing the IceTopSTA3 filter are 
further broken down by the IceTopSTA5 and IceTop_InFill_STA3 filters. The IceTopSTA5 filter requires at least 
5 stations to be hit. The IceTop_InFill_STA3 filter, which aims to select events with energies as low as 100 TeV, 



requires at least 3 stations from the infill tanks (Section 2.2 1 to be hit 



There is an additional IceTop filter, called InlceSMT JceTopCoin (|8J, which is independent of the IceTopSMT 
trigger. It requires the in-ice simple multiplicity trigger to be fulfilled and at least one HLC hit in IceTop. Events 
which pass the InlceSMT JceTopCoin filter consist of coincident events triggering in-ice IceCube with mostly just 
one IceTop station. These events will be used for vetoing high energy inclined cosmic ray showers in addition to 



being used to test and calibrate the entire IceCube detector (see Section 10.6 1 



4.3.5. Data formats and satellite transfer 

For all events passing the filter conditions IceTopSTA3 or InlceSMTJceTopCoin the total charge and the 
leading edge time of all DOM launches is recorded in a condensed format called SuperDST and transmitted via 
satellite. Using a logarithmic representation the IceTop charges, which span 5 orders of magnitude, have to be 
encoded by 14 bit to ensure a precision of about 3 x 10~ 3 throughout the whole dynamic range. A total of 5 byte 
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Table 8: IceTop filters with their inclusive rates. The rates for full waveform transmission via satellite after associated prescaling (4th column) 
and the corresponding data loads are reported in the last 2 columns. All events including the prescaled ones are transmitted in the condensed 
SuperDST format as described in the text. 



filter 


condition 


rate 


full waveform subset 


satellite data 






[Hz] 


prescale 


rate [Hz] 


[GB/day] 


IceTopSTA3 


IceTopSMT & # Stations > 3 


22.3 


10 


2.2 


0.9 


IceTopSTA5 


IceTopSMT & # Stations > 5 


6.1 


1 


6.1 


3.1 


IceTop JnFill_STA3 


IceTopSMT & # In-fill Stations > 3 


3.7 


1 


3.7 


1.3 


InlceSMTJceTopCoin 


InlceSMT & # IceTop HLCs > 1 


100 


100 


1.0 


0.5 



per DOM is necessary for encoding charge, time, station number, DOM number and LC bit. For the observed 
mean number of 22 launches per event the average event size yields 1 10 bytes. The estimated IceTop contribution 
to the data rate for all events in SuperDST format is less than 500 MB per day. 

In addition to the transmission in SuperDST format, the full waveform will be transmitted for selected or 
prescaled events according to the scheme summarized in Table [8] For the IceTopSTA5 and IceTop JnFill_STA3 
filters all events will have the full waveform information, while for the others only some fraction of events will 
have full waveform information for control studies. 



Transient rate monitoring. In 1 s intervals, the rates of the SPE and MPE discriminators (see Section 3.3.3 1 as 
well as of the Local Coincidences are written out for monitoring and detection of transient events. These so-called 
scaler rates are also transmitted for each run via satellite. There is a delay of typically a day until the data are 
accessible in the North. 



5. Calibration 

The tank signals are calibrated in two steps: First, the PMT and the DOM electronics are calibrated in order to 
obtain signal charges in units of photoelectrons and a time reference which is the same for all DOMs in all parts 
of IceCube. This part is common to all DOMs in IceCube. The second step is the calibration of the tank signal 
charges in units "vertical equivalent muon" (VEM) which is the signal a vertical energetic muon generates in a 
tank. 

5.1. Calibration of DOM electronics and PMT 

Calibration of PMT gain and electronic components of a DOM are performed during special calibration runs 
by software called 'DOMCal' running on the DOM CPU (TJ. The basic calibration quantity is the charge recorded 
for single photoelectrons which are produced by cathode noise and radioactive decays in the glass pressure sphere 
(in the following referred to as 'cathode noise'). 

Charge calibration. The charge calibration of the ATWDs starts by converting the ADC counts of each bin to 
voltages. The ADC-count to voltage relation is determined by varying the ATWD reference voltage and applying 
a linear fit to the resulting amplitudes yielding slopes and intercepts for each ATWD bin. The 'pedestal pattern', 
given by the intercepts, is equalized for all bins to a common baseline by the DOM firmware. This common 
baseline is continuously determined during the data taking runs as described in Section [4.2.3| Since the pedestals 
of all ADC cells are equalized, only one number for a channel's baseline has to be obtained by averaging for a 
certain time period, usually for a run of about eight hours duration. The baselines are continuously monitored and 
the calibration constants are updated when needed. 

The ATWD front-end amplifiers are then calibrated in a two-step process. First, the channel with the 16x 
pre-amplifier is calibrated by injecting single photoelectron-like pulses with the pulser on the mainboard, and 
comparing the measured pulse charge to the known true pulse charge from the pulser. The other channels are then 
calibrated consecutively relative to the highest-gain channel using LED light injected into the PMT. 

The PMT gain is calibrated using single photoelectron pulses from cathode noise. The charge distribution 
obtained from these pulses can be described by a Gaussian single-photoelectron peak and an exponential contri- 
bution at low charges 1 1 1 1. The relation between high voltage setting and single-photoelectron peak is determined 
by varying the bias voltage between 1200 and 1900 V. These voltages are higher than those used for the low gain 
PMTs in IceTop (750 V and 1250 V for low gain and high-gain PMTs, respectively, Table [4]i so that an extrapo- 
lation is necessary. The observed single-photoelectron pulses allow the absolute calibration of charges in units of 
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Figure 17: Distribution of muon track lengths in a tank for a cos _ ~ 8 zenith angle distribution. 



PE with an estimated precision of 10%. Note, however, that for IceTop the final precision is given by the VEM 
calibration as described below. 

While for in-ice DOMs the discriminator thresholds are calibrated in terms of single photoelectron charges the 



IceTop discriminator thresholds are set at fixed voltages as described in Section 3.3.3 The voltages can be related 
to a charge in numbers of PE or VEM for average IceTop pulses as demonstrated in Fig. 13 The charge values are 
obtained from the calibrated ATWDs. 

Timing calibration. The pulse transit times from the PMT to the recording are determined using pulses from the 



LED on the DOM mainboard (Fig.[T0|b). 

The master clock of the IceCube experiment is a GPS clock in the IceCube Lab. The local 20 MHz oscillator of 
each DOM is compared to this master clock in a procedure called "reciprocal active pulsing calibration" (RAPcal) 
iflOl at a frequency of 1 Hz. Each signal from any IceCube detector component gets a timestamp with the same 
time reference expressed in UTC times. From experimental data the precision of the time calibration was found to 
be better than 3 ns |Q] . 

The RAPCAL process involves analog measurements on pulses exchanged between the IceCube Lab and 
individual DOMs, and was optimized for the dispersion and attenuation characteristics of in-ice cables. Since the 
IceTop cables are shorter, electronic filters are used to make the corresponding pulse shapes similar to in-ice pulses. 

5.2. Tank calibration using atmospheric muons 

The DOM response to the same energy deposition in a tank varies from tank to tank because of different optical 
properties of the tanks and DOM efficiencies. Therefore, the tank signals are calibrated by normalizing to the signal 



charge of a vertical muon defining the unit 'Vertical Equivalent Muon' (VEM) as described in Section 3.3.1 

The VEM calibration procedure uses the natural flux of atmospheric muons. Since 2009 a single-muon calibra- 
tion trigger is run continuously together with the normal data taking, which has several advantages as compared to 
employing special calibration runs as done before. In particular, there is no detector downtime, and the calibrations 
can be done more regularly. 

5.2.7. Single muon spectra 

The average muon energy at the detector level is in the order of 2 - 2.5 GeV Since muons at these energies 
are minimum ionizing, they always lose about the same amount of energy for a given track length in the tank. For 
instance, the mean energy deposit of a minimum ionizing vertical muon in a tank of 90 cm ice (density 0.92 g/cm 2 ) 
is about 165 MeV Cherenkov light is emitted either by the muon itself or by the secondary electrons along the muon 
track. The total number of Cherenkov photons scales with the energy deposit, and thus with the muon track length, 
in the tank. Therefore, the spectrum of muon signals in a tank is mainly determined by the angular distribution of 



the muons and the tank geometry. Figure 17 shows a purely geometrical calculation of the track length distribution 
in a tank, using a realistic muon zenith angular distribution dN/dD. ~ cos 23 6 Il20ll . The sharp peak at 90cm is 
mainly due to vertical tracks passing through top and bottom of a tank whereas shorter track lengths are entirely 
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due to edge-clippers. Because of the relation between the track length and the number of Cherenkov photons, 



the experimental muon spectrum (Fig. 18 a) exhibits the same features but they are smeared out due to statistical 
fluctuations of the Cherenkov photons, the optical properties of the tanks, and the DOM electronics. Furthermore, 
there is a background under the muon peak mainly from electromagnetic cascade particles from low energy air 
showers. 

5.2.2. Taking calibration data 

The calibration of the IceTop tanks is based on signal spectra from events with preferentially only one muon 
hitting the tank. These spectra are recorded for high-gain DOMs in parallel to the normal data taking by special 
calibration DOM launches, called VEMCal launches, which require the LC condition between the DOMs of the 
same station not to be fulfilled (see section |4.2.3| l. With this anti-coincidence condition, contributions from ex- 
tensive air showers are largely suppressed and those of single muons are enhanced. The DOM resident firmware 
pre-scales the discriminator trigger rate by a factor of 8 192 (= 2 13 ) and launches the DOM readout if there is no LC 
fulfilled, otherwise the following discriminator trigger is taken. This results in a DOM launch rate of about 0.2 Hz 
per DOM. For each of these VEMCal launches the ATWD waveforms are sent to the IceTop DOMHub (Section 



4.3.1 1. When a hub receives such data it generates a calibration trigger which initiates event building around the 
hit (Section |4.3.2| i and passes the data to the online processing. 

VEMCal launches are only enabled for the high-gain DOMs since the gain and discriminator thresholds of 
the low-gain DOMs do not allow triggering on single muons. For low-gain DOMs the VEM calibration has to be 
determined from the overlap of charge measurements by each DOM in a tank, as will be explained below. 

The VEMCal hits are calibrated and extracted during the online processing at South Pole exactly the same way 
as the physics data which are (currently) processed in the North. The raw waveform, given in ATWD channel 
counts, is calibrated using the DOMCal calibration constants (see Section |5~T) i. Corrections for residual baseline 
offset and a droop correction are also applied. Finally, the charge, given in PE units, is calculated by summing 
up all the waveform bins. The extracted VEMCal hits are stored in a special data container and transferred to the 
North within the normal physics data stream, while the original VEMCal events are only stored on tape (unless 
they satisfy other filter conditions). 

5.2.3. Processing of calibration data 

In the North the VEMCal hit data for each individual DOM are collected for one week and analyzed. The 
analysis involves fitting of two plots : The fit to the muon spectra to extract the muon peak of the high-gain DOMs 
and a fit of the correlation between high-gain and low-gain charges to cross-calibrate the low-gain DOMs. 

The muon spectra are fitted with an empirical formula which includes vertical muons, edge clippers and the 
electromagnetic background [21 1: 



/(*) = Po 



. 1.85 1 
L(x;puP2) + 



Pi exp(^f) + l 



+ p3 ■ exp (p 4 ■ x) (7) 

fern 



The first term in the muon part, L (x; p\,pi), describes through-going muons by a normalized Landau distribution 
(from ROOT library 1 22 1), the second the edge-clipping muons by a Fermi-like function and the f em term describes 
the electromagnetic background by a single exponential function with two parameters (p^, pn). The parameter 
Po is the number of muons which are not edge-clipping, p\ and pi are the location and width parameters of the 
Landau distribution, respectively. The number of edge-clippers is 1 .85 po determined from the purely geometrical 
considerations in Fig. [T7J the Fermi function is normalized for positive x values by the factor l/p\. Figure 18 a 
shows a fit of the function |7} to a charge spectrum. 

5.2.4. VEM definition 

The unit 1 VEM is defined as the charge value at 95% of the muon peak which is the maximum of the muon 
contribution in the fit (term in (|7|). The measured muon peak position, obtained with a well-defined fitting 
procedure, is the essential calibration reference. Shifting the VEM definition to 95% of the peak value has been 



chosen because the measurement of tagged vertical muons resulted in a 5% lower muon peak (see Section 5.3 1. 
However, the choice of a 5% shift does not affect physics results because it also appears in simulations. The charge 
1 VEM corresponds to about 125 PE for most high-gain DOMs except for those commissioned in 2005 (the tanks 



are listed in Table [2| which have about 200 PE (Fig. 19 1. The low-gain DOMs have systematically lower VEM 



charges of about 105 PE with correspondingly higher values of the 2005 tanks. This systematic shift between 
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Figure 18: Calibration monitoring plots of the DOMs in tank 61A from April 15, 2012. a) Charge spectrum of the high-gain DOM in Tank 61A. 
The black curve corresponds to a fit employing the function {7) with the muon (green) and electromagnetic (red) contributions. The maximum 
of the muon contribution is indicated by the green arrow while the black arrow points to the value of 1 VEM at 95% below the maximum, 
b) Relative charge difference of the high-gain (HG) and low-gain (LG) DOMs in tank 61 A. The red curve is the average charge difference for 
a given LG charge slice which is fitted by the green curve according to the function jHJ. 
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Figure 19: Histogram of the VEM values of all the high-gain and low-gain DOMs from VEM calibration April 15, 2012. 



high-gain and low-gain DOMs is not fully understood (but accounted for by the VEM calibration), see discussion 
in the next Section l5.2.5l 

The energy of air showers can only be determined by comparison with simulations. It is therefore essential 
that the simulation uses the same definition of VEM to quantify the tank response to air showers. Here the merit 
of the VEM calibration becomes particularly obvious: the details of each tank do not have to be simulated, but one 
typical tank behavior can be used in the simulation, and the simulation can be approximately applied to all tanks 
when signals are expressed in units VEM. The 'calibration' of the simulated signals in units of VEM is described 
in Section |9J] 

5.2.5. Low gain cross-calibration 

Since the muon spectra are only recorded and fitted for the high-gain DOMs the low-gain DOMs are cross- 
calibrated assuming that the two DOMs of a tank record on average proportional light yields per particle. Therefore 
the signal ratio of the two DOMs should be constant in an overlap region where saturation does not play a role. 
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Figure 18 b shows the relative charge differences between high gain and low gain versus the low-gain charge in 
units of PE. The mean values of the scatter plot at a given low-gain charge are constant up to about 2200 PE where 
saturation of the high-gain DOMs starts. 

For each tank the average relative charge differences as a function of the low-gain charge are fitted by the 
expression 

f(x) = (Po for *<P3 (g) 

{ po + p\ log I0 (l + p 2 (x - p 3 )) for x > pi 

with the free parameters po, pu P2 and p-$. The cross-over between the two function parts in |8]l is at x — p^,. 

The horizontal line to the left of the cross-over has an offset po = Ar from zero corresponding to the relative 
difference in the efficiencies of both DOMs which are systematically lower for low-gain DOMs as can be seen in 



Fig. 19 The efficiency includes the quantum efficiency of the PMT as well as the optical efficiencies of the tank 
which may vary within a tank due to inhomogeneities in absorption and scattering. Furthermore a potential mis- 
calibration of the gain of the low-gain DOMs, likely due to extrapolation of the PMT gain below the single-photon 
calibration range (Section [5T| , could be the cause for the systematic shift between low-gain and high-gain DOMs. 
The offset from zero can be converted to a correction factor for the low-gain DOM for the charge of a VEM (or 
any other fixed signal) according to: 

fi LC (VEM) = -i— e„ c (VEM) . (9) 
1 - Ar 

The left side is the VEM value of the low-gain DOM in units of PE if 2hg(VEM) is expressed in units of PE. 
The cross-over at x = pi (in terms of low-gain charge) indicates the beginning of the saturation of the high-gain 
PMT and determines the switching from high-gain charges to low-gain charges for the waveform measurements. 
In terms of high-gain charge (in units PE) the cross-over is at 

QG£" = (i-Mp3- (io) 

5.2.6. Monitoring VEM calibration 

The VEM calibration as shown in Fig. 18 is done for all DOMs in an automated way. Thus stability checks 
and long time DOM response changes can be monitored, as depicted in Fig. [20] for the DOMs in tank 61 A (for 
DOM numbering see Fig. [6]). The variations of the number of photoelectrons per VEM for all DOMs between 
two successive calibrations are within 10% with an rms spread of 3%. The increase of the S /B ratio of the 
muon signal to the electromagnetic background (see definition in Eq. ( fT4] i) indicates increasing absorption of the 



electromagnetic shower component and is related to the increase of the snow height on the tank (Section 6.3 1. 
The discriminator trigger rates plotted below the S /B ratio are mostly determined by the electromagnetic shower 
component and thus decrease with increasing snow. 

5.3. Comparison to runs with muon taggers 

Muon telescope at the South Pole. A portable, solar-powered muon telescope based on scintillation detectors was 
employed in order to tag nearly vertical muons for calibration studies. Comparing the tagged muon events to 
simulations allows verifying the shower generation and the tank simulation. The variation of the tank response as 
a function of muon position and the tank-to-tank variations have also been investigated l23l . 

Measurements were done in the polar season 2005/2006 on tanks deployed one year earlier. The muon tagger 



was placed on top of an IceTop tank, and all hits in both DOMs were recorded for six hours. Figure 21 shows 
the charge spectra for the two DOMs in tank 39B and superimposed the tagged muon charge spectra with the 
muon telescope placed above the center of the tank. The muon telescope spectrum peaks about 5% lower than 
the calibration spectrum because in the latter case the accepted zenith angle range is wider leading to on average 
longer particle trajectories in the tank. These measurements are the justification behind choosing the definition of 



one VEM to be 95% of the muon peak, as described in Section 5.2.4 



Zenith dependence of muon flux. The muon flux at the South Pole was measured for five zenith angles, 0°, 15°, 
35°, 82.13° and 85.15°, with a scintillator muon telescope incorporating an IceTop prototype tank as the absorber 
Il24l . Besides the flux measurements also the zenith angle dependence of the recorded waveforms were analyzed 
and compared to simulations. 

The comparison of waveforms from vertical and nearly horizontal muons showed two features in the data: 
(i) The average amplitude of horizontal muons is about two times as big as that of the vertical muons. (ii) The 
average rise time of the these pulses in the tank is 3 ns longer than that of vertical muon pulses while the decay 
time of the average waveform is nearly the same as that of the average waveform of vertical muons. Both features 
were well reproduced by a Geant4 simulation of the tank response. 
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a) History of DOM 61-61 
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b) History of DOM 61-62 
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Figure 20: Calibration history of the high-gain DOM (a) and low-gain DOM (b) in tank A of station 61 (same as in Fig. |18( . From top to bottom 
are displayed: The number of photoelectrons per VEM, the high-gain - low-gain cross-over used in the calibration of low-gain DOMs, the ratio 
of the muon signal to the electromagnetic background from the fit of the muon spectrum, the discriminator trigger rate (MPE for high-gain 
DOMs and SPE for low-gain DOMs), the high voltage of the PMT and the temperature measured on the DOM board. The green vertical lines 
indicate the time of DOMCal runs where in general DOMCal constants can change. 
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Figure 2 1 : Charge spectra of a high-gain DOM (a) and low-gain DOM (b) in tank 39B collected by the calibration procedure (VEMCal, black 
lines) compared to the spectrum tagged by the Muon Telescope (filled histogram) [23]. 



Laboratory calibration measurements. To study the response of an IceTop tank in more detail a laboratory setup 
has been developed. It comprises an IceTop tank filled with 76 cm of water (not frozen), equipped with two 
optical modules with analog readout. Tiled plastic scintillator plates, mounted around the tank, are operated in 
a coincidence mode. They allow the selection of muon tracks passing through the tank at defined angles and 
positions. For each position and angle, typically 10 4 muons were measured and the mean value and width of the 
distribution of the number of photoelectrons in each optical module was determined. Since the dimensions (76 cm 
height vs. 90 cm) and optical properties of the active medium (water vs. ice) were not the same as at South Pole, 
the aim of the test was (i) to confirm that the results are correctly reproduced by the simulation program 'Tanktop' 



(described in Section 9.2.1 1 which is the reference for the IceTop simulation (Section 9.2 1, and (ii) to gain a better 
understanding of possible systematic uncertainties of the calibration measurement. 

A scan of nearly vertical muons passing through the tank at different x-y-positions has been performed. In 
Fig. [22] a top view of the tiled scintillator centered on the IceTop tank is presented. The left plot shows for each 
segment the mean number of photoelectrons generated by muons which enter through the segment and exit through 
the corresponding segment below the tank. The signal variations due to different distances of the muon track to the 
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Figure 22: Top view of the tiled scintillator centered on the IceTop tank. The mean number of photoelectrons generated by nearly vertical 
muons entering through each segment and measured in an optical module located at (x = -0.285 m, y = 0) is shown. The measurements (left) 
are compared to simulations with the simulation code 'Tanktop' (right). 
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Figure 23: Mean measured number of photoelectrons for muons passing with different track length through a water-filled tank. The track 
lengths are selected under different zenith angles by the scintillator coincidences as described in the text. The same track length can be obtained 
under different zenith angles. The measurements (blue) are compared to simulations. 



optical module, indicated by the small circles, amount to about 40% maximum and are described to within 10% 
by the simulation shown in the right plot. 



Figure 23 shows the mean measured signal strength in units PE for different muon track lengths in the tank. 



The signal strength scales nearly linearly with the effective track length. Since the track lengths have been selected 
by the scintillator coincidences each measurement corresponds to a zenith angular range. The average zenith angle 
for each point is also displayed in the figure. The measurements are well reproduced by the simulation. 

Further investigations of the zenith angle dependence of signal sizes are anticipated, since the test tank has not 
the same geometrical dimensions as the deployed tanks. In addition, a data-based study, e. g. exploiting isolated 
muons on the periphery of well-reconstructed showers, would allow better control of angle-dependent systematic 
uncertainties in physics analyses. 



5.4. SLC charge calibration 



For SLC hits (Section 4.2.2 1 integrated charges and time stamps are delivered by the DOM firmware. Since 



the baseline calibration as described in Section 4.2.3 is not available at this point, the baseline correction is applied 



off-line using the same baseline determination as for HLC hits. Remaining differences to the HLC charges are 
corrected by comparing the correlation between HLC and SLC pulses when both are available. 
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Figure 24: The temporal behavior of the South Pole stratosphere from May 2007 to April 2009 is compared to IceTop DOM counting rate and 
the high energy muon rate in the deep ice. (a) The temperature profiles of the stratosphere at pressure layers from 20 hPa to 100 hPa where the 
first cosmic ray interactions happen, (b) The IceTop DOM counting rate (black -observed, blue -after barometric correction) and the surface 
pressure (orange), (c) The IceCube muon trigger rate and the calculated effective temperature (red) 1 25 1. 



6. Environmental conditions at South Pole 

6.1. Overview 

Since cosmic ray measurements with IceCube use the Earth's atmosphere as a converter medium, variations 
of the atmosphere will affect air shower measurements. The shower development is mainly determined by the 
overburden, X(h) = p{h)/g, at a height h and a pressure p(h) (g is the Earth acceleration, here assumed to be 
constant). Pressure measurements are available at the South Pole from balloon flights. In addition to the total 
overburden, as given by the ground pressure (at South Pole the average is about 680 hPa), also the density profile 
of the atmosphere influences the air showers. A change in density changes the balance between strong interactions 
and weak decays of mesons and thus the rate of lepton production. The high-energy muon rate observed in the 
deep ice is mainly influenced by the density (or temperature for a given pressure) profile in the stratosphere. 

The South Pole atmosphere undergoes a pronounced annual cycle. In the winter months, April to September, 
when the sun is below the horizon, the air is much colder than in summer, with surface temperatures ranging 



from -20 °C down to -70 °C. Therefore, the winter atmosphere is much denser than in summer. Figure 24 shows 
the seasonal variation of the density in the upper atmosphere and the correlation with IceTop DOM rates and in- 
ice muon rates |25|. The high-energy muon rate in the deep ice follows closely these density changes with an 
approximate +(8-9)% variation (Fig. |24p), while the variation of the IceTop DOM rates is about +5% (Fig.[24p). 
In contrast to the strong annual variation of the density profile, there is no cyclic variation of ground pressure. 



Instead changes happen on a much shorter time scale, as seen in Fig. 24 *>. While the high-energy muons rates are 
not influenced by the ground pressure changes, there is a strong correlation between pressure and IceTop DOM 
rates. An increase in air pressure causes a stronger attenuation of air showers leading, for a given shower, to a 
decrease of the observed shower size in IceTop. The related shift of primary energy thresholds makes the IceTop 



rates anti-correlated to air pressure. The pressure effect can empirically be corrected for (blue curve in Fig. 24 b), 
which is useful for rate monitoring. For air shower reconstruction it is necessary to find a correction for the shower 
size S , the observed shower signal at a reference radius, which is the most important observable of an air shower 
(the determination of S is described in Section[8]). 

Another effect which influences the shower development is the snow accumulating on the tanks. While the 
atmosphere affects the overall profile of an entire shower, snow accumulation affects individual tanks depending 
on how deeply they are buried. 
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6.2. Atmosphere 

6.2.1. Barometric effect on rates 

Empirically, the change in the IceTop rates with the surface pressure P as shown in Fig. 24 depends almost 
linearly on the absolute pressure change: 

-^-=BdP => N = N exv(8(P~Po)). (11) 

The barometric coefficients B, introduced by this relation, are regularly determined for the various trigger and filter 
rates by fitting data at different pressures. The fitted B values range from about 0.6 to 0.7 %/hPa for different rates. 
For monitoring purposes a rate N meas measured at a particular pressure P t is normalized to a rate at the reference 
pressure Pq according to 

Nf orr = Nf eas exp(-B(Pi-P Q )). (12) 

After correcting the DOM rates for surface pressure variations (blue curve in Fig. [24 b) seasonal effects are 
clearly visible in the DOM rates, although not as strong as for the in-ice muon rates in Fig. 

6.2.2. Effect of the atmosphere on energy measurements 

At different pressures air showers are differently absorbed in the atmosphere so that an effect on the measured 
shower size is expected. Corrections are derived by comparing measured shower size spectra from periods with 
different surface pressure and finding the correction which makes the spectra coincide. Experimentally it is found 
that the shower size parameter S changes with the surface pressure P with respect to a reference pressure Pq as 

cm crp^ ( p ~ p " \ dlnS(P) 1 .... 

The absorption parameter Ap(S, 9) is studied as a function of the shower size S and the zenith angle 8. For vertical 
showers above 1 PeV one finds l/A p w 0.25 x 10~ 5 Pa~' or, in terms of overburden, l/A p * 0.25 x 10~ 4 cm 2 /g 



The seasonal effect on the shower size measurements was estimated to be A log 10 S ~ ±0.01 for the two extreme 
atmosphere densities in winter and summer ll26ll . Corrections of atmospheric effects on energy determination from 
reconstructed showers are discussed in Section [531 

6.3. Snow on tanks 

While precipitation at South Pole only amounts to about 2 cm per year, there is a significant snow accumulation 
due to drifting. The snow height on top of the tanks increases on average by 20 cm every year with a density of 0.35 
to 0.4g/cm 3 depending on depth. However, accumulation strongly depends on surrounding terrain and buildings. 
At South Pole the wind has a predominant direction and a relatively constant average speed of about 5 m/s.This 
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Figure 26: Snow heights measured in Feb. 2010, Dec. 2010 and Feb. 201 1: a) The ratio S^/Bem versus the measured snow height, b) Compar- 
ison of measured and extrapolated snow heights. The extrapolation is done by using the signal-over-background ratio of the VEM calibration 
with formula jTJ. 



causes an increased accumulation of snow leeward of buildings and slopes. Figure 25 shows the distribution of 
snow on top of IceTop stations in February 2012. 

Snow depths on the tanks can directly be measured only during the Antarctic summer. However, the develop- 
ment of snow heights during the whole year is important for the analyses of the data. For the determination of the 



snow heights between direct measurements the calibration spectra as shown in Fig. 18 are employed. The effect of 



snow is mainly an absorption of the electromagnetic shower component while the muon spectrum remains nearly 
unaffected by growing snow. Therefore, the ratio of the number of muons around the muon peak in Fig. 18 to the 
underlying mainly electromagnetic background depends on the snow height. This ratio is defined by integrating 
over the muon part and the electromagnetic part f em of the calibration fitting function |7]): 

S„ IB EM = Jf """ f M dS I J*'"'" f em dS ( 14) 

Suitable integration boundaries have been found to be S m ,-„ = 0.3 VEM and S max = 2.0 VEM. 



Figure 26 a shows the correlation between measured snow heights and this ratio. The points can be well fitted 
by the function: 

h s = a s x (ln(5^/B em ) + b s ) (15) 

typical values of the parameters are a s w 1.37 m and b s ~ 1.75. 

Formula (JT3J is used to estimate snow heights in between direct measurements. For a correct extrapolation the 
slope a s is fixed to the average of all snow heights but the intercept b s is determined for each tank separately. A 
comparison of snow heights determined in this way with direct measurements is shown in Fig. |26*|b. The fitted 



snow heights are used to account for the snow effects in the reconstruction of shower parameters (Section 8.3 i 
7. Signal processing and data preparation 



Figure 1 1 shows typical waveforms measured in IceTop. For air shower reconstruction, these waveforms are 
processed by a series of software modules to perform three steps: calibration, extraction of total pulse charge and 
time, and pulse cleaning. Only the unsaturated ATWD channel with the highest gain is used in the data analysis. 



First, the ATWD baseline (common for all bins, see Section 5.1 1 is subtracted from the digitized ATWD 
waveforms. The waveforms are then converted from ADC counts into voltage taking into account the gain of the 
ATWD amplifier and the calibration of each individual ATWD bin. In case of SLC hits, where no waveform is 
available, the average ADC calibration of bins 5 to 60 is used. The start time of the waveform is corrected for the 
transit time of the pulse through the PMT and the electronics (Section |5T| i. Finally, all waveforms are corrected for 
the effects of droop caused by the transformer that couples the PMT output to the recording electronics (Section 
3.5 I. The correction uses the inversion of formula ([3]l in Section[33]in a discretized way. This is done by correcting 



the voltage in each time bin for the droop in the same bin and for the summed voltage responses from all preceding 
bins. When a readout window contains consecutive hits from the same DOM, the waveform of the last hit is used 
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Figure 27: Example of the IceTop air shower reconstruction, a) Lateral fit to signals from 25 triggered stations with a reconstructed shower 
size S125 = (65.1 ± 2.8) VEM. b) Time residuals with respect to a plane perpendicular to the shower direction given by Eq. \2l\ . "Upstream" 
and "downstream" refer to tanks being hit before and after the shower core reaches the ground. 



to correct additionally for the residual droop in the follow-on hit. Droop correction cannot be applied to SLC hits 
but is on average accounted for by the SLC calibration (Section |5~4] i. 

The total HLC pulse charge is then extracted by integrating the calibrated waveform and converting it into units 
of PE and VEM to be used in reconstruction. The pulse time is defined by the crossing of the leading edge slope 



of the first pulse with the baseline. The slope is taken between 10 to 90% of the leading edge (see Fig. 11 a). All 
times are expressed in UTC time. 

For each tank, charge and time of the pulse in the high-gain DOM are used, unless the charge surpasses a 



saturation threshold determined during the VEM calibration (expression ( [10] ) in Section 5.2.5 I. In that case, the 
charge is used from the low-gain DOM. The pulse time is always based on the high-gain DOM. If there is no pulse 
in a low-gain DOM within ±40 ns of a saturated high-gain pulse, the pulse is marked as saturated to be treated 
accordingly by reconstruction algorithms. 

In addition, for air shower reconstruction, events are cleaned by requiring the following conditions: 

• A station is discarded if the following condition on the signal times in the two tanks A and B is not met: 

\tA~t B \ < ^ A + 200 ns, (16) 

c 

where and tg are the signal times in the two tanks located at xa and Xg- 

• Stations are grouped in clusters which in principle could belong to the same shower, such that any pair of 
stations i and j in the cluster fulfills the condition 

I*; —Xj\ 

\ti - tj\ < — + 200 ns. (17) 

c 

The station position Xi is the center of the line connecting its two tanks, and f, is the average time of both 
tank signals. In each event, only the largest cluster of stations is kept. 

This selection is done in order to remove obviously unrelated pulses from events. In the IT26 analysis [27| it only 
affected about 4% of events, where on average 2.3 tanks were removed. 



8. Air Shower Reconstruction 

Properties of a primary particle are inferred from the air shower parameters reconstructed from the IceTop 
signals. The reconstructed parameters include the shower core position and direction, and the shower size. The 
latter is a measure of primary energy and is defined as the signal 5 re f measured at a certain distance R K f from the 
shower axis. These properties are reconstructed by fitting the measured charges with a lateral distribution function 
and the signal times with a function describing the geometric shape of the shower front (Fig. 27 1. Currently only 
HLC hits are used in the reconstruction. 
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8.1. Time and charge distribution of air shower signals 

Lateral charge distribution. IceTop tanks are not only sensitive to the number of charged particles, but also detect 
photons via electron-positron pair production. That is different from a scintillator array, which is relatively insensi- 
tive to the photonic part of the signal resulting in a lateral charge distribution which can be described by the NKG 
function [28 , 29 1. Therefore, for IceTop an empirical lateral distribution function was derived from simulations 
MM- 

The charge expectation value in an IceTop tank at distance R from the shower axis is described by: 



R 



-13-K log 10 («/R rer ) 



sm = Sr*-y—\ . (i8) 



which is a second order polynomial in the logarithm of R: 



log 10 S(R) = log 10 5 ref -ft log 10 ( ^- ) - k log? 1 ) . (19) 



^ref) ^V^ief 

This function behaves unphysically at small distances to the shower axis (R < 1 m). However, as described in 
the next subsection, all signals that are closer than a minimal distance to the core, are excluded from the fit. The 
free parameters of the function, in addition to the shower size, S re f, are /3 and k, corresponding to the slope and 
curvature in the logarithmic representation at R — R le f. Simulation studies suggest that k can be fixed without a 
significant impact on the fit result. The current default is k — 0.303. 

At medium energies, the distance from the shower axis corresponding to the average logarithm (log 1() R/m) of 
signals participating in the fit is about 125 m. The current standard fit uses a constant R Te f = 125 m, which was 
found to minimize the correlation between the parameters S Ie [ and f3 on average for the IT26 configuration. The 
shower size parameter is thus referred to as S 125. An optimization of this parameter for the full detector is under 
study. 

Time distribution. The arrival times of the signals map out the shower front. The expected signal time of a tank at 
the position x is thus parametrized as 

t{x) = t + i(* - x c ) ■ n + At(R). (20) 

Here, to is the time the shower core reaches the ground, x c is the position of the shower core on the ground and 
n is the unit vector in the direction of movement of the shower. For fixing the core position, 'ground' is defined 
as the y[S -weighted average of participating tank altitudes, which varies by about 6 m over the whole array. The 
term At(R) describes the shape of the curved shower front as a function of distance R to the shower axis and is the 
time residual with respect to a plane perpendicular to the shower axis which contains x c . A plane through x c is 
described by At(R) = 0, used as start value for the fit, see Section [53] Studies on experimental data showed that 
the shower front can be described by the sum of a parabola and a Gaussian function, both symmetric around the 
shower axis BTl : 

At(R)=aR 2 + bll-exp{-^\\, (21) 

with the constants 

a = 4.823 10~ 4 ns/m 2 , b= 19.41ns, cr = 83.5 m. 

The energy, zenith angle and mass dependence of these parameters is under study and presents an opportunity 
for further improvement of the reconstruction quality. Furthermore, inclined showers will not be symmetric with 
respect to the shower core, which may require a modification of the signal time and charge distribution. 

Function ( [20] > is fitted to the measured signal times with five free parameters: two for the core position, two 
for the shower direction and one for the reference time to- Hence, the complete air shower reconstruction has 
the following parameters: position of the shower core (x c ,y c ), shower direction 8 and 0, shower size S125, slope 
parameter /3, and time at ground to. 



8.2. Shower fluctuations 

Fluctuations of signal size and arrival time for each tank have to be included in the fit procedure. The signal size 
fluctuations are described by a normal distribution of log^S, around the fitted expectation value log 10 Sf l , with 
standard deviations cr logioS depending on the signal charge. The charge dependence of cr logioS has been determined 
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experimentally from the local shower fluctuations between the two tanks of a station and are reasonably well 
reproduced by simulation 1 32 1 . For the shower reconstruction a functional description is used: 



log 10 (cHog 10 s) 



0.283 - 0.078 log I0 (5/VEM) log 10 (5/VEM) < 0.340 

-0.373 - 0.658 log 10 (5/VEM) + 0.158 log 2 lQ (5/VEM) 0.340 < log 10 (5/VEM) < 2.077 
0.0881 2.077 < log 10 (5/VEM) 



(22) 



The arrival time fluctuations are described by a normal distribution with a standard deviation cr,(Rj), 



o-,(Ri) = 2.92 ns + 3.77 ■ 10~ 4 ns • (/?,/m) 2 , (23) 

depending on the distance of tank i to the shower axis as found in experimental data l32l . No energy dependence 
of cr t has been included so far. 

8.3. Likelihood fit 

Likelihood function. The functions ([19}, |20} and ( 121} describing the expectations for the charge and time of air 



shower signals are fitted to the measured data using the maximum likelihood method (Fig. 27 1. In addition to the 
signal charges and times, the likelihood function also takes into account stations that did not trigger so that the full 
log-likelihood function consists of three terms: 

£ = £ q + £o + £ t . (24) 

The first term, 

v (log 10 5 i -log 10 5f t ) 2 v flt 
A = ~2. Yo^ <S*> Z^^^siS^), (25) 

describes the probability of measuring the charges 5, if the fit expectation value at the position of the tank is 5 fit as 
given by the lateral distribution function ( j18} . The standard deviation o"i og|0 s is given by ( |22] i. The sum runs over 
all tanks that have triggered. The second sum in £ q accounts for the proper normalization of the signal likelihood 
and has to be added because the standard deviations depend on the fitted signals. 
The next term of the log-likelihood function d24|, 



£ = 2>(l-(/f) 2 ), (26) 

j 

accounts for all stations j that did not trigger. The probability that one tank in station j delivers a signal at a given 
charge expectation value is 



1 r ( (log 10 S;-log 10 Sf) 



Pf = f exp 

V2W logioS (5f) J 

i°gio s r 



2o%Sf) 



&\2\ 

dlog 10 5,. (27) 



The lower integration limit is defined through the charge threshold 5* r for the tank signal. The charge expectation 
value, 5^', is evaluated for the center of a line joining the centres of the two tanks. Since the two tanks of one station 



operate in coincidence and tanks with unusually missing partners, as described in Section 4.2.4 are removed by 
the cleaning described in Section [7] there are no single untriggered tanks. There is a natural correlation in the 
signal expectation values of two nearby tanks because they have a similar value of the lateral distribution function. 
However, in Eq. ( p6l > the no-hit probability of the are assumed to be independent because the fluctuations about the 
average expectation value, 5^\ of both tanks are uncorrelated. 

The third term of ( f24] >, £,,, describes the probability for the measured set of signal times, 

where the index i runs over all tanks, f, is the measured signal time of tank i, f?' = t(Xi) is the fitted expectation 
value according to function d20|, and cr, is defined in (|2"3"]>. 
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Figure 28: Reconstruction quality of showers in two zenith angle ranges for proton showers with 5 or more stations triggered in the IT73 
configuration. Shown are the 68% resolutions as a function of primary energy, a) Distance between true and reconstructed shower core, b) 
Angle between true and reconstructed shower direction. 



Fit procedure. The likelihood fit is seeded with first-guess calculations for the core and the direction of the shower. 
As a first estimate of the core position the centre-of-gravity of the positions jc, of the m tanks with the highest signals 
weighted with the square root of the charges is calculated: 

xcoc = ■ (29) 

The weight ~\fSj is chosen based on a study of the achievable fit accuracy. The restriction to the m tanks with the 
highest signals (current default is m = 7) yields a better estimate, in particular near the boundaries. The starting 
value for the direction is obtained by fitting to the signal times of the m tanks with the highest signals a plane, that 
means setting At(R) = in the function ( |20] i. 

The likelihood minimization is then done in several iterations to improve the stability of the fit. At first the 
shower direction is fixed and only the lateral fit of the charges is iterated with the free parameters Sn5,P, and core 
position. After each iteration those tanks are removed from the fit that are closer than a minimal distance to 
the shower axis. Iteration is stopped when no more tanks are removed from the fit. The parameter /?„„■„ is optimized 
to improved core resolutions by mitigating the effect of saturated pulses and of biases in the determination of core 
position^] Then, a third iteration step is done in which the shower front curvature is included and the shower 
direction is varied. Finally the direction is fixed again and S125, /3, and core position are allowed to vary for a fine 
tuning. 

The fit performance is demonstrated in Fig.|28]where the energy dependence of the core and angular resolutions 



for the full likelihood fit are displayed. The performance for the energy reconstruction is discussed in Section 10.2 



Small Shower reconstruction. If only 3 or 4 stations are included in the fit, the last two iteration steps are dropped 
because the small number of stations with small separation between the tanks of a station does not offer sufficient 
redundancy. The shower direction and the fo parameter are taken from the first guess fit assuming a plane shower 
front. 

Snow correction in reconstruction. During the fit procedure the absorption of shower particles due to snow accu- 
mulation on the tanks and the snow around the tanks is corrected for. The correction is applied to the fit values in 
each iteration rather than to the measured charges. In this way the 'no-hits' probability for signals due to the thresh- 
olds can be correctly computed using Eq. (|27j». For the correction an exponential absorption model is assumed: 

5r r = £exp(-— ^—V (30) 
\ A s cos 6) 

Here 5, is the charge expectation for the tank i given by the lateral distribution function and S f rr the corresponding 
snow corrected value; h' s is the snow height above tank i, 8 is the zenith angle of the shower axis, and A s is an 
effective absorption length in snow. We are studying if an effective absorption describes correctly the interplay 
between absorption and regeneration in a shower and the different absorption behavior of different particles or if 



'In the IT26 analysis 1271 R m ,„ = 1 1 m was used. 
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Figure 29: Shower size spectrum measured with the IT73 configuration in the zenith angular range cos# > 0.8 in an area of 5.2 X 10 s m 2 . The 
spectrum contains three separate data sets with 3 and 4 stations, 5, 6 and 7 stations and 8 or more stations. 



a more complicated description is required. From studies of data and simulations this absorption length has been 
found to lie between 2 m and 4 m. Dependences of snow attenuation on energy, signal size, zenith angle, and 
position with respect to the core are still under study. 

Atmosphere corrections of shower size. The shower size S125 is corrected for surface pressure dependence by 
inverting formula ( fT3j ). For the analysis of the IT26 energy spectrum 1261 the following corrections were found 
(replacing in ( fT3) l d In S = In 10 log 10 S and the pressure P by the overburden X = Pig): 



0° < 6 < 30° 
30° < 6 < 40° 
40° < 6 < 46° 



dlo 


§10^125 




dX 


dlo 


§10^125 




dX 


dlo 


§10^125 


dX 



= (-1.2 ±0.9)- 1(T 4 



1 



(-9.0 ± 1.5)- 10" 
(-14.3 + 2.9)- 10~ 4 



g/cm 2 ' 
1 



g/cm 2 
1 



(31) 



g/cm- 



The dependence on the atmospheric density profile has not yet been explicitly corrected for. For data of a whole 
year the effect should cancel out. For more restricted time periods different atmospheres are studied in simulations 
(see Section|9~T|. 



8.4. Shower size spectra and energy determination 

A shower size spectrum reconstructed from data taken with the IT73 configuration in 2010/11 is shown in 



Fig. 29 The zenith angular range is cos 8 > 0.8 and the shower core is required to lie within a polygonal fiducial 
area with a size of 5.2 x 10 5 m 2 excluding all border stations. The spectrum contains three separate data sets with 
3 and 4 stations, 5, 6 and 7 stations and 8 or more stations, depicting the different energy ranges. 

The relation of the shower size parameter to the cosmic ray energy depends in a non-trivial way on the primary 
mass composition and the zenith angle and has to be evaluated using simulations (Section [9j. A determination 
of the cosmic ray energy spectrum requires the knowledge of the mass composition. Measurements with IceTop 
alone can exploit the zenith angular dependence of the shower size on the mass Il27l . A better handle on the mass 
composition is achieved by measuring the air shower on the surface, which is dominated by the electromagnetic 
component, in coincidence with the high-energy muons in IceCube. The energy spectrum and mass composition 
have then to be extracted by multi-variate algorithms, for example a Neural Network ll33l l34l or two-dimensional 
unfolding (35). 



To get a rough idea of the corresponding energies of the shower size spectrum in Fig. 29 one may use above 
S125 ~ 1 VEM the approximate relation E/Sns ~ 1 PeV/VEM for proton primaries at small zenith angles. 
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9. Simulation of air showers and the IceTop detector 

The relation between the measured signals and the energy of the primary particle, as well as detection effi- 
ciency and energy resolution of the IceTop detector are obtained from simulations of air showers generated by the 
CORSIKA program [36 37] and then traced through the detector by a Geant4-based Il38l[39l program. 

9.1. Air shower simulation 

Air showers in the atmosphere are simulated using the code CORSIKA. The hadronic component of the air 
showers is simulated using the models SIBYLL2.1 fl40ll4TI or QGSJET-II [42 43 1 for high energy interactions and 
FLUKA 2008.3 fi4ll45l for low energy interactions (the current default is < 80GeV). The electromagnetic com- 
ponent is simulated using EGS4 |46l . 'Thinning' (reduction of the number of traced particles |47|) can optionally 
be applied. 

For the South Pole, CORSIKA provides six different pre-defined atmosphere profiles [37 1. Four of these 
atmospheres (CORSIKA labels 11 to 14) are based on the MSIS-90-E model ll48l . one for each austral season. 
The two others, modeled by R Lipari, are alternatives for austral summer and winter. There is also the option 
to define custom profiles, for example using actual measurements of the South Pole atmosphere by satellites and 
balloons. 



9.2. Detector simulation 

The output of the CORSIKA program which includes the shower particle types, positions and momenta at the 
observation level of 2835 m, are injected into the IceTop detector simulation which uses the Geant4 package. The 
simulation starts with the generation of the amount of light produced by the shower particles in the tanks followed 
by the simulation of the PMT, the DOM electronics and the trigger chain. 

To make the Geant4-based detector simulation more efficient, optical tank properties are studied in a stand- 
alone, single tank simulation (Tanktop). Some of these results are used in a parametrized form for the full detector 
simulation. 



9.2.7. Study of optical tank properties with Tanktop 

A standalone program, Tanktop, has been developed allowing detailed studies of the optical properties of the 
tanks and DOMs. The optical properties are described by wavelength-dependent functions for PMT efficiencies, 
reflectivities of the tank walls and the ice-Perlite or ice-air boundaries and the absorption in the ice. Reflectivities 
of the tank wall liners (Fig. Q and of the perlite filling on top of the ice have been measured under laboratory 
conditions and have to be tuned for the simulation to obtain agreement with observed signal characteristics. The 
tuning was done by scaling the reflectivities by a power law R'(A) = R(A) P , with p typically about 0.38 (that means 
that the tuning requires higher reflectivities than measured in the laboratory). The absorption length in ice was 
taken to be 100 m from measurements in the deep-ice (this value is not critical). For the wavelength dependent 
PMT efficiencies averages of sample measurements are used IfTD . Another study concerned the implementation 
of the snow cover of the tanks into the simulation. 

The simulation studies confirm that an effective, average description of detector properties is sufficient if one 



employs the VEM concept to calibrate data to a standard signal (Section 5.2 1. The essential conclusions from the 
Tanktop studies for the full detector simulation are: 

- an effective exponential decay time of pulses due to reflectivities and absorption; 

- verification that the number of generated Cherenkov photons in the wavelength interval 300 nm to 650 nm is 
proportional to the number of detected photoelectrons; 

- verification that the latter relation is - within the assumptions of the simulation - independent of the direction 
and location of the simulated particle's trajectory in the tank; 

- verification that the effect of snow can be included by starting the shower particle tracking just outside the 
snow volume. 



38 



Table 9: Table of materials and their properties used in the Geant4 simulation. The column 'composition' the fractions are given in parentheses. 



Material 


Composition 


Density [g/cm^] 


Conditions/Properties 


air 


rs{\j. 1 jjzO/ ), 


U.UUlzlO 


temperature T = 243.15 K, 




r\/r\ 11 1 1Q 1 \ 




pressure p — 67 kPa, 




Ar(U.Ulzo2/), 




ionisation potential 








U[ = oj.i ev 


snow 


IT f~\ 


n. to 




ice 




U.Vz 


index of refraction n = 1.31, 








optical photon energies: 








c-opt — [i.yi ev, 4.13 evj 


plastic 


TJ fA 


1 .4zj 


substitute for the tank materials 


Perlite 


N(0.701o), 0(0.2491), 


0.1598 


mixture of 92.92 % air 




ri2<J(U.UU21J, Ar^u.uiiyj, 




and /.Uo 70 reriite 




Si(0.0240),Al(0.0051), 








K(0.0025),Na(0.0024), 








traces of C, Fe, Ca, Mg, Mn 






glass 


S1O2 


2.254 


DOM pressure sphere 


DOM 


SiQ 2 


0.2 


effective DOM material 








with average density 



9.2.2. Tank simulation 

The Cherenkov emission inside the tanks is simulated using Geant4. All structures of the tank, the surrounding 
snow, including individual snow heights on top of each tank, as well as the air above the snow are modeled 
realistically (see Tables [3] and |9]| 49]). The snow heights used in the simulation are implemented for the specific 



data period as explained in Section 6.3 In order to save computing time, Cherenkov photons are not tracked; 
only the number of photons emitted in the wavelength interval 300 nm to 650 nm is recorded. This procedure 
was studied using the stand-alone program Tanktop, described above, which includes Cherenkov photon tracking 
until photons reach the PMT The propagation of Cherenkov photons is modeled by distributing the arrival times 
according to an exponential distribution, which was tuned with Tanktop by varying the optical tank parameters 
such that simulated waveform decay times match those observed in experimental data (Table|6]l. In the simulation 
the number of Cherenkov photons per VEM is fixed and the same for all tanks (Section [93b; the corresponding 
number of photoelectrons per VEM, in turn, is taken from the individual VEM calibration of the real tanks. While 
the photoelectrons are simulated with the exponential time distribution of the pulse decay (Table IS), a realistic 



leading edge slope results from the simulation of transport through the electronics (Section 9.2.3 1. 

To further reduce the required computing time, particles whose trajectories miss a tank's volume by more than 
30 cm are excluded from the detector simulation. Additionally, when the total charge created inside a tank is larger 
than 2000 VEM, further particles are not injected into the Geant4 simulation, but only counted. The final number 
of photoelectrons at the PMT is then scaled accordingly. This significantly reduces the time needed to simulate 
tank signals very close to the core. 

9.2.3. PMT and DOM simulation 

In the next step the generated photoelectrons are injected into a simulation of the PMT followed by the analog 
and digital electronics of the DOM. To simulate the photomultipliers, single photoelectrons with Gaussian shaped 
waveforms (cr = 3.2 ns) and a charge randomly chosen from the measured single photoelectron spectrum IfTTI 
are superimposed. Pulse saturation is simulated depending on the instantaneous current according to the gain 



dependent functions depicted in Fig. 30 (PMT saturation is discussed in Section 3.4 1. 

In the DOM simulation, the pulse shaping due to the analog front-end electronics is applied to the output of 
the PMT simulation. This includes the simulation of the droop effect using a discretized form of formula (3 
(Section 3.5 1 and the individual shaping of the signal paths to the ATWD and the discriminators using formula (T 
(Section 3.3 1. Next, the discriminators are simulated and the local coincidence conditions are evaluated. Finally, 
the digitization of the waveforms by the ATWDs and fADCs and the array trigger (IceTopSMT, Section 4.3.3 1 are 
simulated. 

Simulated data have the same format as the experimental data and are reconstructed in the same way, as 
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1 

10° 10 1 10 2 10 3 10 4 

Figure 30: Saturation curves as used for the IceTop simulation for the high-gain (red) and low-gain (blue) PMTs. The the curves are labeled 
with the corresponding gain value. The saturated current is plotted versus the ideal, unsaturated current. The saturation currents given in Table 
[6]are the unsaturated currents for which the output current is 50% lower. 

described in the previous section. 

9.3. Definition ofVEM units in simulation 

The VEM 'calibration' of the tank simulation, that is the transfer of the experimental VEM definition to the 
simulation, is achieved by relating the number of Cherenkov photons generated by any particle in the tank to a 
value in units VEM. The relation of 1 VEM to the number of photoelectrons is taken from the calibration data of 
a typical tank (for each tank type), see Section [5] However, the specific value of PE per VEM is not important for 
the absolute calibration because it cancels out, but it has an influence on the fluctuations. 

The value 1 VEM is assigned to a certain number of Cherenkov photons such that the shower simulation yields 
a charge spectrum with the correct muon peak position. This assigment is done in an iterative way along the 
following steps: 

- Simulate CORSIKA air showers in an energy range in which showers contribute to the VEMCal spectrum 
with single muons (empirically: 3GeV to lOTeV for an £T 27 spectrum and zenith angles up to 60° |50|) 
and simulate the number of Cherenkov photons emitted in the tank using Geant4. 

- Assign the signal charge of 1 VEM to a certain number of Cherenkov photons (as a start value the Cherenkov 
photons which a vertical muon of some energy produces). 

- Take the relation of PE per VEM from data of a typical tank. 

- Simulate the electronic chain for the photoelectrons. 

- Apply the same calibration trigger condition as in data (VEMCal launch, Section|4]i. 

- Generate the single muon spectra as in the VEM calibration and fit with the same function as in data formula 
{7} in Section|5T2]>. 

- Compare the position of the muon peak in VEM units to the data of a tank (if VEM units are used it can be 
any tank). 

- If the simulated peak is not at 1/0.95 VEM tune the VEM assignment to the number of Cherenkov photons 
(for the factor 0.95 see Section [5T24l >. 

The current implementation of the VEM unit into the simulation is: 1 VEM corresponds to 32 235 Cherenkov 
photons in the wavelength range 300 to 650 nm. This is in Geant4 the yield of a vertical muon with a kinetic 
energy of about 0.8 GeV 
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9.4. Simulation production 
9.4.1. Properties of datasets 

For cosmic ray analysis, output files from CORSIKA simulation are available for showers initiated by different 
primaries (proton, iron, helium, oxygen, and silicon), energies between 10 4 and 10 9 GeV and zenith angles of the 
incoming primaries between and 65°. Simulations up to 40° are used to analyze IceTop events in coincidence 
with the in-ice detector. Above 40°, simulations are employed for studies of surface events and events with core 
locations outside the IceTop array area that can still trigger the in-ice detector (for example veto studies, see Section 



10.5 1. In addition, for energies between 10 7 GeV and 10 9 5 GeV, thinned showers (see below) are simulated. The 



overlap with the energy range of non-thinned showers is used to estimate the systematics induced by thinning. 

The simulation datasets are grouped in energy bins of size 0.1 x log 10 (£'/GeV) with an E~ l spectrum, inde- 
pendently normalized in each bin. For comparison with data, the simulations can be weighted to yield a realistic 
spectrum. The direction zenith and azimuth angles are generated to yield an isotropic distribution (uniformly 
distributed in cos 0, for sampling on a horizontal detector surface the distribution is proportional to cos 6 d cos 8). 

Most of the events are simulated for one specific atmosphere (currently model 12 [37 1 which is the South 
Pole atmosphere of July 1, 1997). Smaller data sets generated with other atmospheres are available for studies of 
atmospheric effects and their corrections. The default interaction models used are SYBILL at higher energies and 
FLUKA in the lower energy range. Additional datasets generated with the QGSJET model instead of SYBILL 
are produced for comparison. Compositions according to various models, such as a two-component model ifSTl 
or the polygonato model [52, 53], are produced by adding different mass components with an appropriate spectral 
weighting. 

9.4.2. Re-Sampling 

Since shower generation is CPU intensive, the same showers are sampled several times before being passed to 
the detector simulation. Showers are sampled inside a circle with an adaptable, energy dependent radius around 
the center of the IceTop array. The number of samples is chosen such that every shower would remain on average 
only once in the final sample after applying selection cuts. This balances effective use of the generated showers 
with the artificial fluctuations introduced by oversampling. Typically showers are re-sampled within radii ranging 
from 600 m to 2900 m for energies above 10 5 GeV (usually no re-sampling at lower energies). 

9.4.3. 'Thinning' 

Thinning of simulated showers [37] is a method to reduce computation time and storage volume by reducing the 
number of particles contributing to the shower development. This is achieved by keeping only one particle out of 
all secondary particles originating from a particle interaction with energies below a threshold energy £,/„„ = exE p , 
where e is the 'thinning level'. The kept particle gets a weight assigned which is the sum of energies of the particles 
below the threshold from which it was selected, divided by the kept particle's energy. This ensures the conservation 
of energy, although not of the particle number. The kept particles can interact again, leading to a multiplication of 
weights. Optionally, to reduce the fluctuations induced by thinning, no more thinning is applied to daughters of 
that particle if a weight exceeds a maximum weight w max . 

For IceTop/IceCube simulations a thinning level of e — 10~ 6 is used for primary energies up to E p = 10 8 4 
GeV. Above this energy the thinning level is calculated as e — 273 GeV /E p so that muons capable of penetrating 
deep into the ice (defined as E > 273 GeV for vertical muons) are not assigned thinning weights. The weights are 
limited to w max - ex E p . 

For the detector simulation 'un-thinning' has to be performed as follows: a weighted particle which hits a 
predefined sampling area around a tank generates a number of clones, each with the same energy as the parent, 
which are uniformly distributed over the sampling area. The number of particles is given by the weight such that 
the energy remains conserved. Technically only those clones are generated which have a chance to hit the tank. 
The radius of the sampling area is optimized such that on average only one particle hits a tank. 

9.4.4. Computing resources 

Cosmic ray simulation is performed on different, geographically dispersed computer clusters which are con- 
trolled by a GRID-based system. The produced data are transferred to the data center at the University of Wiscon- 
sin, Madison and become available for further processing. 

At present, simulated CORSIKA showers for coincident events occupy 145 TB which is a large fraction of 



the available storage. In Table 10 typical CPU times and data sizes for single CORSIKA showers are reported. 
Simulations of the highest energy showers in the range between 10 7 and 10 9 GeV can take up to 2 days for a single 
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Table 10: Some typical CPU times for shower generation with CORSIKA. 



energy 
[GeV] 


not thinned 
time [s] size [MB] 


thinned 
time [s] size [MB] 


10 5 


117 


5 






10 6 


879 


62 






10 7 


8,119 


760 


1950 


167 


io 7 - 9 


60657 


6650 


3314 


260 
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Figure 31: The curves show the average charge in PE per VEM as a function of time between July 2009 and April 2012 for all DOMs which 
have been operational at the given time. The upper curve is for high-gain DOMs and the lower for low-gain DOMs. 



shower with resulting file sizes of up to 10 GB. Simulations in the range from 10 4 to 10 7 GeV result in file sizes of 
the order of hundreds of MB with computing time of the order of hours. 

The detector simulation needs about 1/3 of the corresponding CORSIKA computing time and adds only 3% to 
the disk usage. 



10. Performance 

10.1. Tank properties and detector stability 

The commissioning history of IceTop tanks is reported in Table [2] About seven years after the first tanks were 
commissioned the long-term stability of the tank and ice properties as well as of the electronics can be evaluated. 
In November 2008, two test tanks that had been deployed in November 2003 were inspected to confirm that the 
quality of the ice and the bonding between the DOMs and the ice were still good. More quantitatively the stability 
is demonstrated with Fig. [3T| showing the average signal charge per VEM as a function of time. The high-gain 
DOMs do not exhibit an obvious trend. For the low-gain DOMs a drop of about 10% over the displayed 30 months 
or 4% per year cannot be excluded. The changes per VEM calibration have typically a standard deviation of 
+(4 - 6) PE. The PMT gains are usually adjusted by not more than 15% (the gain precision is estimated to be about 
10%). Only one DOM out of the total of 324 IceTop DOMs shows permanent failure. The uptime of the detector 
is typically larger than 99%. 



10.2. Air shower detection in IceTop 

Effective area. Air showers are reconstructed if 3 or more stations have triggered as shown in Fig. 29 For the 



sub-sample where at least 5 stations have triggered, Fig. 32 shows the effective area as a function of energy from 
simulations of the IceTop array as commissioned in 2010 (IT73). The reconstructed shower cores were required to 
lie within a fiducial area of about 5.2 x 10 5 m 2 with a polygonal shape excluding the outer stations. The effective 
exploitable energy range stretches from about 300 TeV to 1 EeV with nearly 3 x 10 7 events per year observed 
between 300 and 400 TeV, 200 events per year above 300 PeV and about 20 events per year above 1 EeV in the 
zenith range from cos > 0.8. There is no indication of array saturation up to 1 EeV lf54l . Including the 3 and 4 
station events the energy range can be extended to lowest energies of about 100 TeV due to the in-fill stations. 
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Figure 32: Effective area for proton showers which triggered 5 or more stations as a function of energy for the IceTop array as commissioned 
in 2010 (IT73). 
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Figure 33: Energy resolution for proton showers which triggered 5 or more stations as a function of energy for the IceTop array as commissioned 
in 2010 (IT73). 



Energy, direction and position resolutions. The obtained energy, direction and core position resolutions for some 



selected energies are summarized in Table 11 In Fig. 33 the energy resolution as a function of energy is shown 



for proton showers in the IT73 configuration (same data set as for Fig. 32 1. The corresponding plots for the core 



and direction resolutions are displayed in Fig. 28 



These resolutions were determined in simulations by comparing the reconstructed with the generated variables. 
In IceTop, resolutions can also be studied experimentally since the same events can be reconstructed independently 
by the sub-array of A tanks and that of B tanks. From the comparison one can obtain an experimental measure of 
the accuracy of reconstruction. The sub-array resolution analysis is in good agreement with the simulation results. 
However, it should be noted that this method is not sensitive to systematic errors common to the whole detector. 

Signal fluctuations. The sub-array method can also be used to determine the signal fluctuations which have to be 



included in the likelihood fit as described in section 8.3 Signal fluctuations are caused by detection fluctuations 
and fluctuations in the shower development. Both together are measured by comparing signals in the two tanks of 



Table 11: Summary of air shower resolutions obtained with IceTop IT73 for protons in a zenith angle range cost? > 0.82. A resolution is 
defined such that 68% of the events deviate less than the listed value. 
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Figure 34: Separation of elements in the primary cosmic ray composition as obtained with a neural net algorithm. The plot shows for pro- 
ton, iron and a helium-oxygen mixture the separation in terms of a neural network output together with the rms width of the corresponding 
distributions. From [34]. 



a station (with proper corrections for different core distances). The detection fluctuations can be separated from 
the shower fluctuations by running both DOMs of a tank at the same gain and comparing their signals. In [T] it 
was shown that the detection fluctuations are much smaller than the shower fluctuations confirming that the tank 
design provides adequate (or better) resolution. 

10.3. Cosmic ray anisotropy 

Cosmic ray anisotropies have been studied by IceCube with high energy muons in the deep ice Il55ll56l . An 
analysis using IceTop alone is currently in progress l57ll . While the declination angular range is more restricted 
than for IceCube (to less than about 60°) the much higher energy resolution due to the shower reconstruction allows 
finer energy binning than for the muons detected with the in-ice detector. Anisotropy measurements can be done 
in the full IceTop energy range given above; sensitivities to intensity changes of less than 10~ 3 can be reached at 
least up to 10 PeV. The angular resolution of better than 3° is sufficient for studies of anisotropy patterns on a scale 
above about 10°. 

70.4. Primary composition 

10.4.1. Primary composition from coincident events 

A strength of IceCube is the possibility to measure high energy muons in the deep ice in coincidence with the 
shower reconstructed in IceTop. The full IceCube detector can cover the energy range from below the knee to EeV. 
Showers generated by primary cosmic rays in this energy range produce multiple muons with energy sufficient 
to reach the depth of IceCube. For a primary energy of 10 15 eV, for example, proton-induced showers near the 
vertical produce on average about 10 muons with E M > 500 GeV and iron nuclei produce about 20 such muons. 
For higher primary energies, the number of muons increases, and the multiplicity in showers generated by nuclei 
approaches asymptotically a factor of A 0,24 times the muon multiplicity of a proton shower [58 1, yielding about 2.6 
for A = 56. As a consequence of the high altitude of IceTop, showers are observed near maximum so the detector 
has good energy resolution, which is important when the goal is to measure changes in composition as a function 
of energy. 

The first analysis of such coincident data is described in 1341 . The data set is constrained to a subarray of 
the detector and a relatively short time period (about 1 month). An example of the separation of elements in 
this analysis is shown in Fig. [34] The separation increases with energy and is in the order of the width of the 
distributions. The measurement of the average logarithmic mass, (In A), from such distributions as presented in 
11341 is still dominated by systematic uncertainties, mainly from energy reconstruction in IceCube from changing 
environmental condition in IceTop and the hadronic interaction models for the shower development. In future 
analyses we expect to reduce the systematic uncertainties and to improve the separability by refining the algorithms 
and including additional mass sensitive variables. Mass sensitive observables are, for example, muon counts at 
the surface (see below), shower shape variables or in-ice observables like stochastic energy loss. A combined 
evaluation of such variables will allow constraining models. 
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Figure 35: Muon counting in IceTop: Distribution of tank signals for various cuts on the signal expectation S exp in the energy range between 1 
and 30 PeV. The average number of muons per event can be statistically extracted by fitting the muon signal together with the background (the 
example shows a fit for S exp < 0.125 VEM). 
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Figure 36: Limits on the ratio of diffuse gamma ray and cosmic ray fluxes from within 10° of the Galactic Plane. The magenta line is the 
90% confidence level limit obtained from one year of data taken with the IC40/IT40 configuration. The limits are compared to an analysis 
by CASA-MIA experiment at lower energies [61] and to the projected sensitivity of 5 years data taking with the full IceCube detector. The 
projected sensitivity is given for the full covered energy range by the dashed blue line and for smaller energy bins by the blue points. 



10.4.2. Primary composition: Muon counting 

Muonic and electromagnetic signals in IceTop can in general not be distinguished. However, muons show up 
as a relatively constant signal of about 1 VEM at larger distances from the shower axis where the expectation value 
of a tank signal, S exp , becomes small compared to a muon signal (S exp is obtained from the fit to the lateral shower 



distribution). As shown in Fig. 35 the muon signal becomes more and more prominent when requiring smaller S exp 
ll59l . For S exp < 0.125 VEM the figure illustrates that the number of muons can be well fitted. Comparing these 
muon numbers as a function of energy to simulations of different primary masses, an independent information on 
the mass composition is obtained. 

10.4.3. Primary photons 

IceCube can efficiently distinguish PeV gamma rays from the background of cosmic rays by exploiting in- 
ice signals coincident with an IceTop event as veto against hadronic showers. Gamma-ray air showers have a 
much lower muon content than cosmic ray air showers of the same energy. Candidate events are selected from 
those showers that lack a signal from a muon bundle in the deep ice. Results of one year of data, taken in the 
2008/2009 season when the detector consisted of 40 strings and 40 surface stations (IC40/IT40) are presented in 
l60l . Figure 36 shows the IC40/IT40 limits on diffuse gamma ray fluxes together with a projected gamma-ray 



sensitivity of the final detector (IC86/IT81). The limits are compared to a CASA-MIA analysis at lower energies 
iffjTl . Point source fluxes have also been searched for and excluded at a level of about 10~ 18 - 10~ 17 cm _2 s _1 TeV _I 
depending on the location in the sky. The complete IceCube detector will be sensitive to a flux of more than an order 
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time relative to the shower front in ns SLC multiplicity 

Figure 37: IceTop veto against cosmic ray background in the search for extremely high-energetic neutrinos studied using data from one month 
running in the IC79/IT73 configuration, a) The plot shows the time distribution of any SLC hit in IceTop relative to the shower front expected 
if the event was a cosmic ray event (see details in text). The different distributions are for in-ice tracks which would pass the nearest IceTop 
tank at the distances given in the legend. The vertical solid black lines define the on-time window (cosmic ray signal), the dashed red lines the 
off-time window (background from random IceTop hits), b) The plot shows for the > 500 m curve in a) the distribution of the number of hits in 
the on-time window (black) and the off-time window (red). 



of magnitude lower, allowing the search for PeV extensions of known TeV gamma-ray emitters. IceCube/IceTop 
has currently the highest gamma ray sensitivity in the energy range between about 1 and 10 PeV. 

70.5. IceTop activity as veto for IceCube events 

The production of extremely high-energetic (EHE) neutrinos is expected in collisions of cosmic rays with 
the cosmic microwave background. The detection of these cosmogenic neutrinos with IceCube |62| requires the 
ability to discriminate very rare and energetic signal events from an abundant background of cosmic ray induced 
muons. High energy cosmic ray air showers produce high numbers of muons densely packed around the shower 
core trajectory. These bundles of muons emit large amounts of Cherenkov light in the IceCube detection volume. 

In 1 63 1 a study is presented which explores the possibilities of using the IceTop array to improve background 
rejection while keeping a large fraction of the neutrino signal. In a limited zenith range around the vertical this 
is obviously possible with high efficiency. However, by exploiting the IceTop SLC hits (single tank hits without 



additional conditions, see Section 4.2.2 1 air showers can also be efficiently rejected at larger zenith angles where 
the shower core does not hit IceTop. 



The plots in Fig. 37 include data corresponding to 33.4 days from runs distributed over the whole year 2010 
(IC79/IT73 configuration). Events were selected as neutrino candidates if they had more than 10 4 photons detected 
in IceCube, corresponding to a neutrino energy of more than 10 17 eV. Under the assumption that the candidate is a 
cosmic ray event, a hypothetical shower front can be constructed from the in-ice track and its movement over the 
IceTop array can be followed. The plot in Fig.[37]a shows the time distribution of any SLC hit in IceTop relative to 
the expected shower front. The different distributions are for in-ice tracks which would pass the nearest IceTop tank 
at different distances of the shower core to the nearest IceTop tank. With the condition that the track passes through 
the in-ice detector the distances of the shower axes from the nearest IceTop tank correspond to zenith angle ranges: 
250m, 500 m and 750m correspond to approximately 10-30°, 30-50° and 25-60°, respectively. The vertical solid 
black lines define the on-time window containing the cosmic ray signal, the dashed red lines the off-time window 
which yields the background of random hits in IceTop. The plot in Fig. [37]b shows for the > 500 m curve in a) the 
distribution of the number of hits in the on-time window (black) and the off-time window (red). 

In the off-time window the hit multiplicity distribution drops very quickly, while the distribution of coincident 
hits is much flatter. The off-time hit distribution can be well described by a Poisson distribution with a mean of 
about 0.25//is. That means requiring for a neutrino event to have less than 1 or 2 SLC hits yields neutrino signal 
reductions by 22% and 2.6%, respectively. The veto efficiency depends strongly on the required amount of light 



detected in IceCube and on the the distance of the extrapolated track to the nearest IceTop tank. From Fig. 37 b, 
where more than 10 4 photons in IceCube and a closest approach of more than 500 m are required, one can read off 
a veto efficiency of about 75% for less than 2 SLC hits. Requiring 10 5 photons in IceCube, the same criteria result 
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Figure 38: The difference between the zenith angle defined by the line connecting triggered IceTop station with the the center-of-gravity of 
triggered in-ice DOMs and the track reconstructed by the in-ice detector. In IceTop it was required that only one station has triggered and that 
the signal is made up by more than 400 photoelectrons. 



in a veto efficiency well above 98%. Exploiting the IceTop veto the cuts on neutrino candidates can be widened 
gaining some 10% overall acceptance. 

10.6. Calibration oflceCube with IceTop 

Events reconstructed with IceTop that are also seen in the deep detector can be used to verify event reconstruc- 
tion and the consistency of the relative alignment of both detector components. Examples of verification of timing 
and direction with IceTop are given in Q. 

Showers that trigger IceTop produce usually bundles of several (at 1 PeV) or many muons in the deep detector. 
In contrast, much of the atmospheric muon background in deep IceCube consists of single muons. Most of these 
events are from cosmic -rays with primary energies below lOTeV. By selecting a sample of coincident events in 
which both tanks at one and only one IceTop station are hit, it is possible to discriminate against high-energy events 
and find a sample enriched in single muons. Such a sample can be selected by the filter InlceSMTJceTopCoin 
(Table[8]in Section 4.3.4 1. Coincidences involving only an interior IceTop station provide a sample in which about 



75% are single muons in the deep detector 1641 . 

An analysis of such coincident events using data taken in 2006 with the IC9/IT16 configuration is presented in 
11641 . Figure [38] shows the distribution of the differences between the zenith angle defined by the line from the hit 
station to the center-of-gravity of hits in the deep detector and the direction obtained from the muon reconstruction 
algorithm in the deep detector. The estimated angular accuracy of the connecting line is about 3°. The observed 
rms of the zenith difference distribution of 3.4° means that the in-ice reconstruction of muons has an accuracy of 
better than 2°. The mean of the distribution being close to confirms the accuracy of positioning of both detectors 
as quoted in Section |2~7o] 

70.7. Heliospheric physics with IceTop 

The IceTop tanks detect secondary particles from multi-GeV cosmic rays with a counting rate exceeding 1 kHz 
per detector. The detector response to secondary cosmic rays is discussed in Ifl8l . 

A certain class of transient events, like solar flares or gamma ray bursts, can be detected with very good time 
resolution via a common rate increase in all IceTop tanks. The observation of the Dec 13, 2006 Sun flare event 
11651 demonstrates these abilities of IceTop. Following this observation the detector readout has been setup such 
that counting rates could be obtained at different thresholds (Section |3.3.3[ ) allowing the unfolding of cosmic ray 



spectra during a flare (see Fig. 14 1. 



11. Summary and Outlook 

The IceTop air shower array, the surface component of the IceCube Neutrino Observatory, is operating in its 
final configuration since the beginning of 2011. In all aspects it has reached design performance. With an in-fill 
array the energy threshold is even lower than in the original design, leading to an energy coverage from about 
100 TeV to 1 EeV. The detector is primarily designed to study the mass composition of primary cosmic rays by 
exploiting the correlation between the shower energy measured in IceTop and the energy deposited by muons in 
the deep ice. 
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First results on the energy spectrum and mass composition have been obtained with the still incomplete detector. 
We look forward to composition results from the complete detector. These results should lead to a substantial 
improvement of our understanding of cosmic rays in an energy range where the transition from galactic to extra- 
galactic origin is expected. 

In addition, the IceTop array has found various other applications. Counting rates of individual tanks as a 
function of different thresholds allow investigation of transient events from the sun and maybe other astrophysical 
objects. Using single tank hits IceTop can be very efficiently used for vetoing cosmic rays in the search for high- 
energy neutrinos. IceTop is also employed for direction and time calibration of the in-ice detector. 

For the future an extension with radio antennas is being studied to enhance the air shower detection capabil- 
ities ||66| . With the measurement of the shower development in air a radio detector would add complementary 
observables sensitive to the composition of cosmic rays. 
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